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We review the recent literature on modified theories of gravity in the Palatini approach. 
After discussing the motivations that lead to consider alternatives to Einstein's theory 
^ ' and to treat the metric and the connection as independent objects, we review several 

' topics that have been recently studied within this framework. In particular, we provide 

, an in-depth analysis of the cosmic speedup problem, laboratory and solar systems tests, 

QQ ■ the structure of stellar objects, the Cauchy problem, and bouncing cosmologies. We also 

' discuss the importance of going beyond the f{R) models to capture other phenomeno- 

logical aspects related with dark matter/energy and quantum gravity. 
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1. Introduction 



5h 

' Einstein's theory of general relativity (GR) represents one of the most impressive 

exercises of human intellect. It implies a huge conceptual jump with respect to 
Newtonian gravity. The idea of gravitation as a force acting in an absolute space 
is replaced by a geometrical theory of space and time in which the space-time it- 
self is a dynamical entity in interaction with the particles and fields living in it. 
This interaction is prescribed by a minimal coupling of those fields to the space- 
time metric according to what is today known as the Einstein equivalence principle 
(EEP). The dynamical equations for the gravitational field itself were deduced on 
grounds of mathematical simplicity and demanding that certain conservation laws 
were satisfied. Unlike the currently established standard model of elementary par- 
ticles, no experiments were carried out to probe the structure of the theory. In 
spite of that, to date the theory has successfully passed all precision experimental 
tests. Its predictions are in agreement with experiments in scales that range from 
millimeters to astronomical units, scales in which weak and strong field phenomena 
can be observecil. The theory is so successful in those regimes and scales that it 
is generally accepted that it should also work at larger and shorter scales, and at 
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weaker and stronger regimes. 

This extrapolation is, however, forcing us today to draw a picture of the uni- 
verse that is not yet supported by other independent observations. For instance, to 
explain the rotation curves of spiral galaxies, we must accept the existence of vast 
amounts of unseen matter surrounding those galaxies. Additionally, to explain the 
luminosity-distance relation of distant type la supernovae and some properties of 
the distribution of matter and radiation at large scales, we must accept the existence 
of yet another source of energy with repulsive gravitational propertieJ^. Together 
those unseen (or dark) sources of matter and energy are found to make up to 96% 
of the total energy of the observable universe! This huge discrepancy between the 
gravitationally estimated amounts of matter and energy and the direct measure- 
ments via electromagnetic radiation motivates the search for alternative theories of 
gravity which can account for the large scale dynamics and structure without the 
need for dark matter and/or dark energy. In this sense, we note that the Newtonian 
accelerations felt by stars and gas clouds in the disk of spiral galaxies are orders 
of magnitude smaller than the accelerations measurable in laboratory. Thus there 
is no experimental evidence supporting the validity of Newton's law down to such 
tiny scales. For this reason, it seems legitimate and well justified to explore modifi- 
cations of Newton's law and Einstein's theory to see if they can provide a consistent 
alternative picture of the observed Universe. In this directions we find Milgrom's 
proposal of Modified Newtonian dynamicsS' (MOND), which fits surprisingly well 
many observational data in its natural regime of applicability and may be related 
with the Palatini approacl3, and a variety of fully relativistic theories which can be 
used to do cosmology, such as f{R) theoriei^, scalar-tensor theories, scalar-tensor- 
vector theoriei^, higher-dimensional and brane-world scenarioJ^, . . . 

The extrapolation of the dynamics of GR to the very strong field regime indicates 
that the Universe began at a singularity and that the death of a sufficiently massive 
star unavoidably leads to the formation of a black hole singularity. Space-time sin- 
gularities signal the breakdown of the theory, because the absence of a wel l-defined 
geometry implies the absence of physical laws and lack of predictabilitji^l^J Por this 
reason, it is generally accepted that the dynamics of GR must be changed at some 
point to avoid these problems. A widespread belief is that at sufficiently high ener- 
gies the gravitational field must exhibit quantum properties that alter the dynamics 
and prevent the formation of singularities. In this sense, a perturbative approach 
to quantum gravity indicates that the Einstein-Hilbert Lagrangian must b e sup- 
plemented by quadratic curvature terms to render the theory renormalizabl 4^ ^ ^ ^ ^ i 
More recent approaches to quantum gravity, such as string theory, also regard GR 
as the low energy limit of a t heo ry that should pick up increasing corrective terms 
at higher and higher energieJ^. The canonical quantization of GR using the so- 
called Ashtekar-Barbero variableJ^ predicts that the continuum space-time of GR 
is replaced by a quantum geometry in which areas and volumes are quantized in 
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bits of an elementary unit of order the Planck scale. The low energy limit of this 
theory should also recover the classical dynamics of GR with corrections signaling 
the discreteness of the space-time. 

The above discussion shows that there are theoretical and phenomenological 
reasons to explore the dynamics of alternative theories of gravity. Though dark 
matter and dark energy could play in cosmology a role similar to that played by the 
neutrino in the process of radioactive beta decajQ, we must try to figure out if our 
theory of gravity can be suitably corrected to explain the dynamics at large scales. 
Since we have well grounded reasons to believe that gravity must be modified in the 
ultraviolet regime, we should not be surprised by having to add corrections also at 
some infrared scale. 

1.1. The Palatini approach to modified gravity 

Because there are no limits to imagination, one should use experiments as a guide to 
constrain the range of possibilities to build an alternative theory of gravity. In this 
sense, the ex per imental efforts carried out in the 1960's to understand the nature 
of gravitation^^ and the kind and properties of the fields associated to gravity, left 
it clear that gravitation is a geometric phenomenon. This led to the conclusion that 
the matter and the other non-gravitational fields must couple only to the metric, 
which implies that the total action must be of the form 

S = SG[gal3, <p,Af„...]+ Srn[9al3, i^] , (1) 

where gai3,4>, A^, . . ., represent the gravitational fields (which can be scalars, vec- 
tors, and tensors of different ranks) and -0 represents collectively the matter fields. 
This defines a class of theories known as "metric theories of gravitjH^' which, by 
construction, should satisfy the EEP. 

Perhaps motivated by the restrictions imposed by the EEP, alternative theo- 
ries of gravity have traditionally focused mainly on pseudo-Riemannian geometry, 
thus forcing the affine connection to be metric com pati ble. However, metricity and 
affinity are a priori logically independent conceptJ^^ and, therefore, there is no 
fundamental theoretical reason to constrain the connection to be metric compati- 
ble. In fact, since affinities are very simple and fundamental geometrical entitieil^, 
in applying Ockham's razor to the construction of alternative theories of gravity 
we should give them higher priority than to other types of tensorial fields. For 
this reason, in this work we will mainly focus on modified theories of gravity in 
which metric and connection are regarded as independent fields. Note that the ex- 
istence of a metric in the theory naturally endows space-time with a Riemannian 



^Since conservation of energy and momentum was a pillar of special relativity, rather than propos- 
ing a modification of this principle, Pauli postulated the existence of a massless particle, the 
neutrino, to explain the spectrum of energies in the process of beta decay. Though in that case a 
"dark matter" particle solved the problem, to explain the anomalous perihelion shift of Mercury 
Einstein had to modify Newton's theory of gravity. 
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connection, the Levi-Civita connection L'^y Thus, accepting the existence of an 
additional connection T^^ is equivalent to having an independent rank-three tensor 
field B'^^ = - L'^^ in the action. 

Though gravitational redshift experiments do not impose very tight constraints on 
the possible coupling of a non-metric connection to matter-, for simplicity we will 
follow the guide provided by experiments and will assume that freely falling bodies 
do follow geodesies of the metric. Thus, rather than working in a purely metric-affine 
framework in which matter is allowed to couple to the independent connection, we 
will consider only this restricted version, which is known as Palatini^ formalism. 
This way we stick ourselves to the class of metric theories of gravity introduced 
above, in which the matter action is only coupled to the metric (and perhaps to 
its derivatives via the Levi-Civita connection) and the gravitational sector is of the 
form SG[gai3,Tp^,<P,A^,...]. 

The Palatini method to obtain the field equations of GR was introduced by Ein- 
stein himself in 1925 . Despite considering independent variations of the metric 
and the connection, the resulting equations in GR turn out to be equivalent to those 
obtained doing variations of the metric only (metric variational formalism). This is 
so because the equation for the connection simply establishes its compatibility with 
the metric. However, this is just an accident. For other Lagrangians, in general, the 
field equations in metric and Palatini formalisms are different, as we will see in de- 
tail later (see Refs ll9l [?T] for some studies on the relation between Palatini and 
metric formalisms). But the differences between metric and Palatini formalisms go 
beyond the field equations, and this can already be seen in the context of GR. In fact, 
since the Einstein-Hilbert action contains second-order derivatives of the metric, to 
have a well defined variational principle one must add a surface term proportional 
to the extrinsic curvature, which explicitly refers to an embedding of the space-time 
into some background metric. In the Einstein-Palatini action, however, there are no 
derivatives of the metric and we only find first-order derivatives of the connection. 
As a result, it is usually claimed that no surface terms are necessary. However, 
to have a consistent formulation, conserved Hamiltonians at infinity in asymptoti- 
cally flat spacetimes, and to correctly reproduce the thermodynamical properties of 
black holes, it has been recently found that a certain surface term must be added to 
the action. This surface term does not refer to any background, b ut whe n there is 
a background available, the metric and Palatini descriptions matctP^I^. This im- 
plies that the corresponding path integral formulations of these two theories may be 
quite different. It is also worth noting that the consideration of the Einstein-Palatini 
action instead of the Einstein-Hilbert one was crucial for the implementation of the 
non-perturbative canonical quantization of the theory using Ashtekar variableJ^. 
Therefore, the Palatini approach must be seriously considered not only to explore 
new phenomenological extensions of GR aimed at explaining the large scale struc- 



''For a discussion of this terminology see Ref llSI 
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ture of the universe, but also as a potential way to make contact with quantum 
gravity phenomenology. 

1.2. Goal and structure. 

In this article we will review the recent literature on modified theories of gravity 
framed within the Palatini formalism. Most of it will deal with theories of the f{R) 
type, in which the gravity Lagrangian is given by a function of the scalar curvature 
R, but we will also comment on scalar-tensor theories and extensions of the f{R) 
family that include other curvature invariants such as R^,yR^'^ . Though the Pala- 
tini approach had been considered in the past in different contexts, the interest in 
Palatini f{R) theories was boosted by the observatiorPSl that some of the known 
problems of the modePSl f[R) = R — /i proposed in metric formalism to ex- 
plain the cosmic speedup, could be avoided by considering its Palatini version. Since 
then numerous works have addressed different aspects of Palatini theories including 
the late-time and early-time cosmology, solar system and laboratory tests, stellar 
structure, the Cauchy problem, black hole thermodynamics, . . . We have tried to 
provide a comprehensive and careful review of the literature on those topics that 
have received more attention. However, since we may have missed some useful and 
important references, we encourage the reader to help us in this task letting us 
know about those works. By "careful" we mean that we have tried to be as precise 
as possible in stating who did what and when, avoiding diffuse lists of references, 
which should help the newcomer find its way through the growing literature. We 
acknowledge that this procedure has been difficult in many cases and also that our 
criterion may not have been the best one at some points. Comments and suggestions 
in this respect will also be very welcome. 

The content has been organized as follows. In Secl2]we provide a detailed deriva- 
tion of the field equations of Palatini f{R) theories and discuss their scalar-tensor 
representation. Since space is limited, for the field equations of other Palatini theo- 
ries we refer to the corresponding literature. Then we split the chronological evolu- 
tion of the literature in subjects: the cosmic speedup problem is reviewed in SeclSl 
laboratory and solar system tests are discussed in Sec|H some questions related 
with stellar structure and new results regarding the Cauchy problem are analyzed 
in SecO and the relation of Palatini theories with quantum gravity phenomenology 
is discussed in SeclHl We end with a summary and future perspectives. 

2. Field equations for Palatini theories 

Since most of the recent literature on Palatini theories has focused on f{R) theo- 
ries, we present here a detailed derivation of the field equations for this case. For 
extensions to actions containing other curvature invariants and couplings to scalar 
fields, we will refer to the corresponding literature. In this section we also comment 
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on the scalar-tensor representation of f{R) theories, the conservation of energy and 
momentum of matter, and some literature on other Palatini theories. 



2.1. f{R) theories 

The action of Palatini f{R) theories is as follows 



(2) 



where Sm is the matter action, represents collectively the matter fields, is a 
constant with suitable dimensions (if f{R) — R, then = SttG), R = g^^R^iv, 
R^f = R^npu, and R^jSfiu — df^T^n — d„T 



^U^tp - r^.r;;^ represents the 



components of the Riemann tensor, the field strength of the connection PJ^^- 



Note 



that since the connection is determined dynamically, we cannot assume any a priori 
symmetry in its lower indices. This means that in the variation of the action to 
obtain the field equations we must bear in mind that PS ^ P"o. We will assume 



a symmetric metric tensor g^^, = g^,^ (for theories with non-symmetric see, for 
instance, Ref. [27| . The variation of the action ([2]) with respect to the metric and 
the connection can be expressed as 



SS = 



1 
2"^ 



xy/-g 



f 



Sg^" + iRg^'^iRpu 



ss„ 



(3) 



where /r = df /Or, and R[^iv) represents the symmetric part of Rpu. Straightfor- 
ward manipulations show that can be written as 



liu - A vp 

X — pA pA 

connection. The contribution of the 8R 
the following expression 



Ap) 



(4) 



where 2^;, = P^, 



P^p represents the torsion tensor, the antisymmetric part of the 



term, I — J d'^xy'^^fRgf^''SRf^„, leads to 



A 



-gJ 



x\ 



srtp {-Va (V^fRg'-') + Vp {V^fRr") + 2V^/flff^-5L}] , 

(5) 



where = }r {g'^'ST^^ - 5^^(5P^^). Having in mind that^ V^V^? = 5^%/^ - 
V^^^Z—g, we find that 

Va (V^jJ^) = 9a (V^J^) + V~9fR [ff^^^^A - 5'ig^''S^,;\ ^P^^ . (6) 

Inserting this result in ([5]) and assuming that the surface term J d^xdx {\/^gJ'^) 
vanishes at the boundaries, the field equations can finally be written as follows 



fRR{p.u) - 7;9piy = f^^T^ 



(7) 



-Va (V^/fl5^'')+^5;Vp {V^fR9'''')+2V^fR {g^^S^^ - S^^'S^^ + g^^'^S'^,) = H 

(8) 



UfJ, 

A 



January 21, 2011 1:5 WSPC/INSTRUCTION FILE PA2MG'V1-1 



Palatini Approach to Modified Gravity 7 

where T^^ = — and H'^'^ = —6Sm/ST^^ = because we assume that the 
matter is not coupled to the connection. To proceed further, it is common in the 
literature to impose the torsionless condition S^^ = 0, which eventually turns ([5]) 
into the simpler form 

Va {V^fRg^") = . (9) 

However, with a bit of extra effort, we will gain deeper insight into the role and 
properties of the torsion and will see that an expression analogous to (O can be 
reached without imposing any restriction on the torsion tensor. The first step is to 
trace over v and A in ([8]) to get 3Vp {y/^fng'^P) = A^/^fng^'PS^p. We then insert 
this result in (|8]) to obtain 

- V;, iV^hr'') + 2V^fR (^g^'^S^,^ - ^-fg^^S^p + g'^'^^^) = (10) 

We now split the connection into its symmetric and antisymmetric parts, which we 
denote C^^ and 5^,^ respectively, and reexpress Va (\/— ff/itff^'') in the form 

Va {^fB.gn = V? {V^fngn+^fn [^''"^L + g^^S^ + g'^'^S^,] , (11) 
where {\/~g f bSJ^^) only depends on the symmetric part of the connection, which 
means that V^A^ — d\Af^ — C^^Ap. Inserting this result in (ITU|) . we get 

V? {V^fngn = V^fR (g^'^^L - g'^'Sl + g'^'^Sl^ - ^-Jlg^^^J,) . (12) 

Adding and subtracting to this equation the same expression but changing the order 
of /i and u we find the following relations 

V? {V~gfRgn - V~gfR [g'^s^^ ~ \ i^W + Kf) si^ (i3) 
- g^'^sl^ = \ {8W - s^yp) s^p . (14) 

Written in this wajQ, it is clear that the symmetric part of the connection is coupled 
to the antisymmetric part (the torsion) via the contraction S'^p. This term is also 
sourcing the right hand side of the torsion equation ([TH) . This fact suggests a new 
step aimed at simplifying the structure of ((T5| and . Consider the new variables 

f ^ = + aStS'^^p , (15) 

and take the parameter a = 2/3, which implies that = T^^^^j is such that 

— 0. The symmetric and antisymmetric parts of the connection F^^^ are related 
to those of F^^ by 

1 



= Cl. + ^ {^IS^p + 6ls:^) (16) 
1 



'^Note that Eqs. ^ and ^ represent sets of 64 independent relations which are equivalent to the 
40 relations of ifTsI l plus the 24 relations of |(T4)| . 
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Rewriting ([T5)) and using these new variables, we find 

Vf {V^fRg'^'') = (18) 
g'-^S^xa - 9"^ SI ^ . (19) 

Written in this form (fT9|) iniphes that Spxa = Saxp, where gjivS^a = Sp\a- Since 
the torsion is antisymmetric in the last two indices, the symmetry of the first and 
third indices automatically implies that S^aq — -i^ S'^^ ~ 0. Using this result in 
(IT7)) we find that 

= I {S^^S^a. - 5^3^,^) . (20) 

This result indicates that the torsion is generated by a vector yl^ = S^^ , which has 
important consequences and will be useful to solve ([T8| . The fact that 

r^. = C^. - 1^,5: (21) 



implies that R$^,{r) ^ R'^^,{C) - |c)[pA,](5;^, from which we get i?^,(r) = 
-^/^Qi/lr) = Rfii^iC) — From this it follows that the symmetric part of 

the Ricci tensor that appears in Q is insensitive to the torsion because (F) — 
i?(^y')(C). Obviously, R is also insensitive to the existence of this type of torsion, 
i.e., i?(r) = R(C). This property is known as the projective invariance of the scalar 
curvature^^ We are now ready to solve for [and ©]• Though seems 

to involve up to second-order derivatives of the connection, it can be reinterpreted 
using the trace of (O, 

RfH - 2/ = K^T . (22) 

This equation implies that R — -R(F) = R{C) can be solved algebraically in terms 
of T, thus leading to i? = R{T) and Jr = fR{T), which are functions of the matter 
and possibly of the metric but not of the independent connection. The solution of 
(|18p can thus be easily found by defining a new metric hfj^i, = fR{T)gfiu in terms 
of which that equation becomes simply Vj(v'— ft./i^'') = 0, which is an algebraic 
equation linear in the connection that leads td^H 

C^u = id^hp, + d,hp^ - dph^,) . (23) 

This completes our analysis of the equations that determine the connection V^^^. 
We have found that, in general, F^^ is made out of a symmetric part, C^^^ plus 
a vector-like contribution — |(5;J^v4^. This vector is responsible for the existence of 
torsion, S*^^ — | {^'^^u — S^A^^^, but it does not affect the metric field equations 
([7|) [see also Eq.(|25|) below], which justifies the usual approach in the literature of 
setting it to zero from the beginning. From this analysis it follows that the four 
conditions = 5*^^ = are enough to force the total vanishing of the torsion. 
When matter is coupled to the connection, the constraint S'^p — has also been 
suggested as a way to avoid potential inconsistencies of the field equations due to 
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the projective invariance of the scalar curvature in Gr[^ and in f{R) theorieJ^Ol 
(see also RefH 



A non-trivial choice for the torsion vector can be motivated by introducing the 
expression ([^ for C^^, in (PT|) . Using the relation /i^^ = fng^v one finds 



r^. = L%+y^ Kd^R ~ <?M-a"/fi] - -S: [A^ - jj^d^Hj , (24) 

where we have denoted = {df_igpi, + d^gp^ — dpg^^). Since the dynamics is 
insensitive to the presence of the vector Ap = Ap — -^^dpfji, one may wish to 
set Ap = to simplify the form of (l24t . By doing this, one finds that S^pi, — 
^{S^S^ - S^d^) and r"^, = L'^p, + K'^p,, where K'^p, = S'^p, + Sp'^^ + S,'' p 
is the so-called contorsion tensor. One can easily check [using for instance Eq. (jlO|) ] 
that this connection turns out to be compatible with the metric gpu, i.e., it ver- 
ifies {y/^gg^'^) ~ 0. This result shows that a torsionless f{R) Palati ni theory 
is dynamically equivalent to a metric-compatible f{R) theory with torsiorP^EHSSl^ 
At the same time, those two particular cases are dynamically equivalent to a non 
metric-compatible Palatini f{R) theory with arbitrary torsion generated by a vector 
field, which is the general case discussed here. 



dpfndyfR ~ \^gpv{d!Rf 



Now that we have an expression, Eq. ()2ip . for the connection in terms of the 
metric, the matter, and the vector A^, we can insert this solution for in ([7]) to 
obtain an equation that only involves the metric gp^ and the matter: 

2 JR. i}R ^UR) 

^ [V,,V,/fl - gp.DfR] (25) 

JR 

where Rp^{g), R{g), and ^ p^ufn are computed in terms of the Levi-Civita con- 
nection of the metric g^i,, whereas R and //j must be seen as functions of T. To 
make our notation clearer, since hpn and gpn are conformally related, it follows that 
R = R{T) = gf^-'Rp^ir) and R{g) = g'^'Rpuig) are related by 

R = ^(5) + TT^^A/i^aVfl - (26) 

^JR JR. 

where, recall, Jr — fR.{T) is a function of T . It is important to note that in vac- 
uum, Tpi, — 0, the solution of ([2^ is just a constani|f| Ryac = RiO), which implies 
that /h(0) is also a constant. As a consequence, the derivative terms on the right 
hand side of (j25l) vanish and that equation boils down to Gpi, = —Aeffgp^, where 



Equation i22l could have more than o ne solution, which could be interpreted as corresponding 
to different realizations of the Universe ESI. For simplicity, we assume that there exists only one 
physical solution, though one should bear in mind that particular models could have various 
solutions that were in agreement with observations in a certain regime. 
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Aeff = ^(.'l ^ plays the role of an effective cosmological constant. This 

means that the dynamics of Palatini f{R) theories departs from that of GR with a 
cosmological constant only in regions that contain sources, where is no longer 

constant and the dfn terms are not zero. Therefore, it naturally follows that outside 
of the sources the solutions take the same form as those of GR with a cosmological 
constant, Birkhoff's theorem holdJ^^, there are only two propagating degrees of 
freedon|2£l, and there are no instabilities'^^ of the kind found in the metric version 



of these theorieJ^Sl. Note, however, that the conditions that such solutions must 
satisfy at the boundary separating the sources from the vacuum region will not, in 
general, be the same as in GR because the interior dynamics is different. 

For some purposes, it may be useful to express the Palatini field equations 
psp using the auxiliary metric h^i, instead of g^^,. Taking into account that R = 
g''''Rf,,,{r) is related with R{h) = h>''' R^,^{T) by R{h) = R/ fn, can be put as 



G^. {h) = ^^^T,. - k{T)h^, , (27) 

where A(T) = {RJr ~ /)/2/| = (/ + K^T)/2p^. It is worth noting that the confor- 
mal transformation needed to go from the representation (^5]) (the so-called Jordan 
frame) to the representation ([TT)) (Einstein frame) has absorbed all the terms with 
derivatives that appeared on the right hand side of (j25p . which makes simpler the 
manipulations of the field equations (|27l) . If one decides to forget about the original 
physical motivations that led to construct the /(i?) theory in the Jordan frame 
and chooses to interpret the Einstein frame metric h^^i, as the physical metric that 
defines free particle geodesies (which implies a redefinition of physical observables) 
then (|27p can be seen as a theory with a density-dependent effective Newton's con- 
stant and a varying cosmological constant A(T). This possibility has also received 
some attention in the literature'^-'^ 42 , 43 , 44 , 45 |^ 



A final comment regarding the vacuum equations is in order. It is easy to see that 
the connection (|23|) is invariant under a constant rescaling of the metric hap — Xhap 
and that G^i^{haf}) — Gf^^{\hajj). If we now compare equations ([25]) and \21\ in 
vacuum, we find that G^^{g) = -Aeffg^i, = -Ae///ip^, with A^ff = /fl(0)Ae//. 
For the discussion of local experiments and stellar structure, it turns out to be 
convenient to rescale the metric in such a way that = h^^, in vacuum. This 
is simply achieved by taking = /"(j) hp^v This leads to Ae// = A^fj, which 
simply states that both constants are measured in the same units. This simple 
observation makes it clear that the difference in the dynamics of Palatini f{R) 
theories in Einstein and Jordan frames amounts to a matter-induced local rescaling 
of units, i.e., the units used in Einstein and Jordan frames differ by a factor that 
depends on the local energy-momentum density. With the constant rescaling, g^^^, = 
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j^^h^i, — j^hfi^, the field equation ([27| can be put as follows 



MT)h 



where = ^7/^(0), and A(T) = {.f/fn.{0) + k^T)/2(P^. 



(28) 



2.2. Scalar-tensor representation of f{R) theories 

The equations of motion (|25p derived above can be rewritten as those of a usual 
(metric-compatible and torsionless) Brans-Dicke scalar-tensor theory, 



1 

2^ 



c^R{g)- {d^<jyd>^^)-V{<P) 



+ S'm[g^i.,V'mI?9) 



by just introducing the following notational change 

= , V{c^) = Rfn - f (30) 

where in order to express V = V{(j)) we assume invcrtiblcEl the relation — fj^ to 
obtain R — R{(p). The equations of motion psp for the metric then become 



1 



(31) 

where in our case the constant parameter w takes the value uj = —3/2. In the 
Brans-Dicke theory, the scalar field </) is governed by the following equation 



(3 + 2uj)a(j) + 2VU) - (j)^ = K^T 

d(p 

which using w = —3/2 boils down to 

2V{(j)) - (j)^ = K^T . 



(32) 



(33) 



This equation is the same as (|22l) but written using the notational change in- 
troduced in ([501) . It is interesting to note that f{R) theories in metric formal- 
ly 

ism also admit a Brans-Dicke-like representation in which w turns out to be 
w = 0. To our knowledge, the identification of Palatini f{R) theories with the case 
w ~ —3/2 was fir st car ried out in Ref. |48l though a scalar-tensor representation 
was already 

knowiPEOl. 

The extension of this result to the metric-compatible f{R) 
case with torsion was first given in Refi33[ and to the more general case discussed 
here in Ref.34, (sec also Refsi35l [511 [52[ for related works). Though this scalar- 
tensor represei itat ion can be useful for some considerations like the computatiorJSSl 
and discussion'^ of the Newtonian and post-Newtonian limits and black hole 



®Note that, unlike other derivations of the scalar-tensor representation, our manipulations do not 
impose any constraint on fan. See, for instance, Ref. 1461 for further details on this. 
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thermodynamic it should not be taken beyond its natural context. In fact, though 
one may be tempted to interpret Palatini /(i?) as the limiting case w —3/2 of the 
Brans-Dicke theorjl^, the fact is that the theory corresponds exactly to the precise 
value w = —3/2. For this reason, the absence of dynamics for the corresponding 
scalar field (absence of the □</) term) is not an issue of fine tunning, and the relation 
between it and the matter needs not necessarily be interpreted as a strong coupling 
regime in which matter is forced to satisfy certain constraints to avoid exciting the 
□0 term. If the scalar field equation is read in its original f{R) form, its meaning 
and implications are much more transparent. Equation (I22p means that geometrical 
objects such as the scalar R are algebraically related with the matter sources in a 
way that depends on the form of the Lagrangian f{R). In GR this relation is linear, 
R = —K^T, but in other theories it may be non-linear, R = R{T). As we will see, 
that relation may end up imposing constraints on the geometry, which obviously 
may back-react conditioning the dynamics of the matter fields. This interpreta- 
tion is naturally extended tomore general Palatini theories which do not admit a 
scalar-tensor representatiorl^SET] Therefore, in the Palatini version of f{R) theo- 
ries, unlike in the metric formalism, the independent connection does not introduce 
new dynamical degrees of freedom. Rather, it modifies the way matter generates the 
space-time curvature associated with the metric by generating new matter terms 
on the right hand side of the field equations. 



2.3. Conservation of energy and momentum 

In Palatini f{R) theories, like in all metric theories of gravity of the form ([T]), the 
conservation of the energy-momentum tensor is naturally sa tisfied and follows from 
the invariance under diffconiorphisms of the matter actiortSSEHSOl, This can be 
seen as follows. Consider the variation of the action under an infinitesimal change 
of coordinates Sx^ — e^(a;) 

2 J Sg^y 

Since the (canonical) energy momentum tensor is defined as T^" = -^^^ — — , 
and a diffeomorphism induces a change in the metric of the form dg^n, = 2V(^e,y), 
where is the usual derivative operator involving the Christoffel symbols of the 
metric g^j/, it follows that 

SS,n ^ \ I d^xV^T'^'^V^e. . (35) 



5S„. = ^ I . (34) 



2 _ 

If the matter action is invariant under diffeomorphisms, SSm = 0, then an integra- 
tion by parts leads to 



SS„ 



i y d^xV^V^T'^"^) = . (36) 
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Since (v^^^'") = y^gV^T'"'' and SS^ vanishes for arbitrary e^, dM]) im- 
plies that V^T'^'^ — 0. Note that despite this elementary result, it has sometimes 
been claimed that Palatini f{R) theories do not satisfy the conservation of energy- 
momentunJ^. The covariant conservation of energy-momentum in very general 
modified theories of gravity was studied in detail in Refi62t in which theories with 
non-minimal couplingJ^Sl between the matter Lagrangian and the curvature were 
also considered. 



2.4. Other Palatini theories. 

In the introduction we emphasized the fundamental role that the connection should 
play in the construction of alternative theories of gravity. In this sense, it is remark- 
able that the consideration in the recent literature of f{R) theories in metric for- 
malism was naturally followed by their Palatini counterpart. However, even though 
they are equally justified, scalar-tensor and higher-dimensional theories (to name a 
few) in the Palatini approach have not received the same attention, and the litera- 
ture in these subjects is scarce. We just find some studies of the conditions in which 
metric and Palatini scalar-tensor theories lead to the same field equationi^SI^ on 
how dimensional reduction in 5— dimensional Kaluza-Klein theory compares with 
the metric approaclf^, an attempt to unify gravitation and electromagnetism in a 
5— dimensional quadratic curvatur e m odel^, and some applications to inflationary 
cosmology^^ and its perturbationJ^ in f{R) and scalar-tensor models. 

The Palatini approach has been recently used by Milgrom^ in the context of 
modified Newtonian dynamics (MOND), which could open new avenues for the 
phenomenology of Palatini theories in the context of dark matter. Milgrom's recent 
approach consists on expressing the Lagrangian formulation of Newtonian gravity 
using a Palatini approach and then introducing the necessary modifications to im- 
plement the MOND equations. This is done considering in the Lagrangian density 



^^ = ^ - 20Vg) +p(-v^-A (37) 



where p — '^^miS{x — Xi), independent variations with respect to the variables g 
and (j). Variation over g yields g = —'S7(f), and over (p gives Vg = —k'^p/2, which 
yields Newtonian dynamics. A MOND-like theory is obtained by introducing the 
acceleration scale aq ^ 10^^° m/s^ and replacing g by a^Q{g^ / a^) in ([57)1 . where 
the function Q{z) must be such that Q{z) — )• z+constant for large z, to recover the 
standard Newtonian laws at high accelerations, and Q{z) ~ (4/3) z'^/^ for z <C 1, to 
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produce MOND at low accelerations. The resulting equations are 



(38) 



(39) 




(40) 



where v{z) = dQ{z)/dz. The above equations indicate that particles move according 
to the standard Newtonian law of inertia in the potential (p. However, as it follows 
from ([301), the MOND acceleration fieli^ ja/ owd = — V0 turns out to be an algebraic 
function of the Newtonian acceleration field g. It is worth noting, as pointed out in 
ReflU that one can get the gravitational part of the Lagrangian ((37|) from the non- 
relativistic limit of the Palatini formulation of GR. In this sense, we want to remark 
that since connections play the role of gravitational accelerations, as is clearly seen 
from the geodesic equation c?it^/o?T + r^^M"u^ = 0, whereas the metric is related to 
the potential field via goo ~ — 1 + 2(/), the failure to satisfy the standard Newtonian 
acceleration law could be seen as a manifestation of connection-related effects. For 
generalizations and relativistic extensions of the theory presented here, see Ref|4] 
(see also section 16.21 for some related results) . 

3. Cosmic speedup in Palatini f{R) theories 

Observations of the cosmic microwave background (CMB) radiation^^, high redshift 
supernovae surveyJ^, large scale structurJ^, and baryon acoustic oscillationi^ 
suggest that the expansion history of the universe has passed through a number of 
phases, which consist on an earlier stage of rapidly accelerated expansion (known 
as inflation) followed by two periods of decelerated expansion dominated by the 
presence of radiation and dust (matter without pressure), respectively, and a current 
phase of accelerated expansion that started some five billion years ago following the 
era of matter domination. The field equations of GR in a Friedmann-Robertson- 
Walker (FRW) spacetime with line element ds^ ~ —dP + a^dx^ filled with non- 
interacting perfect fluids of density pi and pressure Pi, 



where K is the spatial curvature, p = J^iPi^ ^^'^ P — J^i^i^ indicate that a 
phase of positive accelerated expansion can only happen if there exists some mat- 
ter/energy source that dominates over the others and whose equation of state satis- 
fies Px/px < —1/3, where Px and px represent the pressure and energy density of 
that source. A natural candidate to explain the current phase of cosmic acceleration 
is a cosmological constant A, for which Pa/ Pa — —1- However, this simple proposal 




a 



a 




(41) 



Strictly speaking, to recover MOND one should impose the further constraint g = 



-V<t>N- 
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is hard to accept from a theoretical point of view. If A represents a new fundamental 
constant of Nature, one could expect new physical phenomena at cosmic scales in 
analogy with what happened when the Planck constant was discovered. If it is seen 
as vacuum quantum energy, then it is generally claimed that its observed value is 
too small to be in agreement with a naive estimation from quantum field theory, 
though if we apply more rigorous techniques of quantum field renormalization in 
curved space-times the predicted value turns out to be much smalleJ^ than the 
observed one. For these and other reasons, there seems to be a widespread desire to 
explain the current cosmic speedup by means of some dynamical entity rather than 
by a pure constant of cosmic nature. 

The fact that the field equations of Palatini f{R) theories in vacuum exactly boil 
down to those of GR with an effective cosmological constant turned these theories 
into a very natural candidate!! to explain the cosmic speedup. For suitable choices 
of the function f{R), it could happen that the new gravitationally-induced matter 
terms that appear on the right hand side of p5|) were negligible during earlier 
phases of the expansion history but became dominant at later times, thus allowing 
an expansion that closely resembles GR in the past but produces cosmic speedup 
today. One could thus explain the transition from a matter dominated universe 
to an asymptotically de Sitter accelerated one with standard sources of matter 
and radiation but without the theoretical problems posed by a strictly constant A. 
The most famous f{R) model of this kind investigated in the Palatini approach 
was borrowed from a proposal of Carroll et al. ■ in metric formalism, namely, 
f{R) = i? — /i^/i?, where = fi^ /k^ represents the energy-density scale at which 
the effects of the modified dynamics are relevant. VollicH^ considered this model 
and showed that after the standard matter-dominated era, the expansion approaches 
a de Sitter phase exponentially fast. To see this, consider the modified Friedmann 
equation corresponding to a given f{R) Lagrangian in a universe filled with matter 
and radiation 



where pm represents the energy density of the (pressureless) matter, pr is the energy 
density of radiation, and i? is a function of only because T — —pm- In the 1/R 



sit should be noted that the Palatini dynamics is radically different from that corresponding to 
f{R) theories in metric formalism. In that case, the modified dynamics is due to the existence 
of an additional effective scalar degree of freedom which is non-minimally coupled to the scalar 
curvature. This co uplin g turns the metric version of f{R) theories into a particular type of extended 
quintessence model^^land, therefore, the metric f{R) predictions are indistinguishable from that 
type of dark energy models. 




(42) 
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model, one finds 

which recovers R w n'^pm when Pi_i/pm 'C 1 and tends to the constant value Rvac = 
\fZp? when p^/pm ^ 1 (see FigH]). We thus see that when p^/ Pm ^ 1 then 
dHJ behaves as ~ -f^ci? ~ >^^{Pm + ipr/3){p^/pm)'^ + . . ., which is virtually 
indistinguishable from GR. However, when the matter energy density, p,„ ^ 
drops below the constant value p^, p^/pm » 1, then (|^ goes like iJ^ « + 
+ . . ., which tends to a constant and implies an asymptotically de Sitter 
expansion, thus confirming the late time cosmic speedup (see Fig[T]). 




Figure 1. Comparison of the time evolution of the curvature (left) and the expansion factor (right) 
in GR and the 1/R model with the initial condition a{t)\^_-^ = 1. We took p^/pmo = 0.7/0.3. In 
the l/R theory the curvature tends to a constant at late times, thus implying a de Sitter phase. 



3.1. Cosmological constraints 

The 1/R model was soon compared with observations of type la supernovai^, 
though such first studies, as we will see, were excessively optimistic about its vi- 
ability. This optimism may have its origin in earlier studies of Palatini f(R) cos- 
mologies which concluded that these theories were very poorly constrained^, being 
|/flfl(0)| < 10^^'^ one of the constraints coming from cosmological data. Besides the 
R— p'^/R theory, which represented a small departure from GR at low matter densi- 
ties, some authors also explored whether radical departures from the GR dynamics 
at cosmic scales such as f{R) = or f{R) = alni? could be compatible with 
observations. These models were confronted with the Hubble diagram of typ e la 

Supernovae, the data on the gas mass fraction in relaxed galaxy clusters^, and 

77 

baryon acoustic oscillations ' ' . Though the fits to the data were good, the statisti- 
cal analysis did not suggest any improvement with respect to the standard ACDM 
model. On the other hand, tight constraints on the family of models R — aR^ were 
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obtained by studying the cosmic microwave background (CMB) shift parameter and 
the hnear evohition of inhomogeneitiej^ plus the Hubble diagram of type la super- 
novae and baryon oscillationa^. Besides finding that the f{R) = R — i-t^/R model 
was strongly disfavored by the data, it was found that the combined observational 
data were capable of reducing the allowed parameter space of the exponent /3 to 
an interval of order ~ 3 x 10"^ around fi — 0, with a having a value similar to 
the cosmological constant. This meant that R — aR^ w i? — a — a(3 In R could be 
restricted to a tiny region around the ACDM model. More stringent constraints on 
this model were found comparing its predictions with the CMB and matter power 



spectraPSI, pushing the (3 parameter to the range ^ 10~^, thus making this model 
virtually indistinguishable from ACDM. These conclusions have been reconfirmed 
by considering updated datct^ ^'^ and strong lensing statistici^HSSI, A different 
class of models^, with f{R) = (i?" — Rq)^^", has also been confronted recently 
with various data samples. The constraints on the parameters, n = 0.98 ± 0.08, also 
place this model in the vicinity of the ACDM model. 

The models considered so far modify the gravitational dynamics at late times, which 
turns out to be strongly constrained by observations. Modifications at early times 
should be very weak because of the strong constraints imposed by big bang nucle- 
osynthesis and CMB primary anisotropies. One could thus consider whether mod- 
ifications at intermediate times could be in agreement with observations. A model 
proposed in this directiorP3 takes the form f{R) — R + AiiJQe^l^'/''^^^''^ , where 
Hq represents the current value of the Hubble parameter, Ai measures the magni- 
tude of the departure from GR, and A2 controls the time at which the correction 
becomes relevant. Note that at late times this f{R) recovers the ACDM model 
(which corresponds to the limits i? — >■ or A2 — J> 00). Though the background evo- 
lution of this model is not significantly different from the standard ACDM model 
for A2 — 500, 1000, which means that it can hardly be constrained by type la super- 
novae data, its effects on the CMB and matter power spectra are dramatic, being 
A2 = 1000 safely excluded. The strongest constraints are imposed by the matter 
power spectrum. This can be unders tood by lo oking at the growth equation for the 
comovmg energy density fluctuationi^HIMllSl for large momentum k 

where x — loga(t), and = /fl,/(3/fl,(2/fli7 + f^)) represents the effective sound 
speed squared. If > 0, the perturbations oscillate instead of growing, whereas for 
cl < they become unstable and blow up (this happens for f{R) = R — aR^ if 
/3 > 0). In the ACDM model — 0. The form of the matter power spectrum in 
the exponential and power-law models, therefore, changes significantly with time 
developing an intricate oscillatory structure for larger k that clearly confiicts with 
observations, which allows to strongly constrain the parameter space of the models. 
The most optimistic constraints restrict the parameter A2 to the regiorP^ > 
5 X lO'^. 
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In parallel to the c onside rations of above, a theoretical consistency check using 
phase space analysiiSnHII -^^as also carried out to determine whether some families 
of f{R) models could allow for the different phases in the expansion history of the 
universe suggested by observations. It was shown that radiation, matter, and de 
Sitter points exist irrespective of the form of the function f{R) provided that the 
function 

r(Ti\- ^^^^ ~ '^DRIb.r f.^. 

''^''^ '\Rfn~f)iRfnn-fR) ^''^ 

does not show discontinuous or divergent behaviors. Thus models satisfying the 
condition C{R) > —3 lead to a background evolution comprising the sequence of 
radiation, matter and de-Sitter epochs. From this it follows that, unlike in met- 
ric formalism, theories of the type f{R) = R — /3/i?" do allow for the sequence 
of radiation-dominated, matter-dominated, and de Sitter eras if n > —1. For the- 
ories of the type f{R) — R + aR"^ — one finds that an early inflationary 
epoch is not followed by a standard radiation-dominated era, which confiicts with 
the idea that early and late time cosmic acceleration could be unified with this 
type of modelJ^. In particular, for m > 2, the inflationary era is stable and pro- 
hibits the end of inflation; if 3/2 < m < 3, then inflation ends with a transition 
to a matter-dominated phase, which is then followed by late time acceleration; for 
4/3 < m < 3/2, inflation is not possible; and for < to < 4/3 one can have the se- 
quence of radiation-dominated, matter-dominated, and late-time de Sitter without 
early-time inflation. 



4. Solar system and laboratory tests 

Most of the f{R) models found in the literature have been proposed to address 
phenomenological issues related with the largest scales. It is generally argued that 
at such scales the theory of gravity could depart from GR, implying that GR should 
be seen as an approximation valid only at certain scales. This, in a sense, justifles the 
study of f{R) theories that are very far away from GR, i.e., that are not of the form 
f{R) = i?+small corrections. However, this viewpoint should be supported by an 
explicit mechanism able to explain why/how the gravity Lagrangian should/could 
change its form depending on the scales involved. Since such a mechanism has not 
been seriously discussed in the literature on Palatini theories, we assume that the 
proposed f{R) models should be treated in the same way as GR and, therefore, to 
be viable they should agree with observations and experiments on all scales. For 
this reason, the same f{R) models that have been proposed to explain the cosmic 
speedup should be in accord with the dynamics of the solar system and laboratory 
systems. In this section we address these points in detail. 
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4.1. Solar system 

In section [21 we remarked that the field equations of Palatini /(i?) theories in 
vacuum boil down exactly to those of GR with a cosmological constant. For this 
reason, if one considers a non-rotating, spherically symmetric star like the sun, the 
metric outside the star can be written as a Schwarzschild-de Sitter solution 

gf.^dx'^dx'' = -A{r)dt^ + + r^dn^ (46) 



Air) 



with A(r) = 1 — 2GMq/i — Ae//r^/3, where Ag// represents a cosmological constant, 
G is Newton's constant, and Mq is identified with the mass of the star. The well- 
known model f{R) — R — n'^/R, like any other f{R) model, admits such solutions. 
In this particular case, it is easy to see that Ag// — •\/3/i^/4. One is then tempted 
to conclude that this model is compatible with solar system observationJ^ because 
for sufficiently small Ag/ / its predictions are virtually indistinguishable from those 
of the SchwarzschilcPSi and KerJ^ solutions of GR, which pass all observational 
tests. However, the situation is more subtle because, due to the modified dynamics 
within the sources, the transition from the interior solution to the exterior solution 
is not, in general, as simple as in GR. To illustrate this point, let us consider a 
presureless body such as a rocky planet or a gold sphere, for example. For such 
objects a fo rmal analytical solution for an arbitrary Lagrangian /(i?) can be easily 
obtained^ by writing the field equations in the form ([25)1 and taking the ansatz 
ds"^ = gf^^dx'^dx'^ = 



'^h^^dx^dx'^ with 



1 



-B{r)e^'^^-'^de 



1 



where we have defined (/)(T') = 
(see the discussion above Eq, 



to guarantee that outside the sources 



B{r) 

/fi(0) 

1^). We then find 

2 d$ _ ^2 /T, 
r dr <j)^ \ 
d{r[l -B]) _ k^T* 



dr' 



(47) 



= h 



n 



B 



dr 



-A(T) 



(48) 
(49) 



Defining now B{r) = 1 — 2GM{r)/r in (jlH]), we can rewrite M(r) and $(r) as 



/ dx 

2G Jo 


x' 


\n 

02 




(50) 


— dx x 
2 Jo 


± J. 


b^B 




(51) 



If we consider a point outside of the sources at radius r, the above equations can 
be readily integrated leading to 



M{r) = Mo + 
^ ' 6G 

$(r) = $0 , 



(52) 
(53) 



January 21, 2011 1:5 WSPC/INSTRUCTION FILE PA2MG'V1-1 



20 Gonzalo J. Olmo 

where A(0) = /(0)/2/i^(0), Mq, and <J>o are constants. Since we are assuming a 
presureless fluid, = —p, taking units such that = SttG we find that Mq and 
$0 are given by 

dr^^^^l r^drKi-p) (54) 

$0^(5/'%. (55) 

io 02(^„2GM(r)) 

where i?© is the radius of the object (where p vanishes). With these results, outside 
of the sources (where (/)(0) = 1) the line element (|47p coincides with pS)) if we ab- 
sorb the constant factor e^*" in a redefinition of the time coordinate, identify A{r) 
with B{r), and take Ae// = A(0). 

Once an /(i?) Lagrangian is specified, eqs. (|47)) . ([50]) . and (ISTI) provide a complete 
exact solution for a presureless, nonrotating, spherical object. The usual GR ex- 
pressions are recovered by just taking cj) — \ and A = 0. In particular, one finds 
the Newtonian expression for the mass can thus use this 

general solution to study the Newtonian limit corresponding to such objects. We 
start by writing down the general expression for the gu component of the metric 

,„.-^^[,-^^yi'^-^-'.^ (56, 

If we consider, for instance, the f{R) — R — p^ / R theory, using (j43l) we can write 



m = 1 - 



2 

1 + 12 ( ^ 



(57) 



From this expression we see that 4>{T) varies continuously from (p^o = 3/4 in- 
side matter (p ^ p^) to i^o = 1 in vacuum. This should have disastrous conse- 
quences for the theory because in the solar system the gu component of the metric 
is only slightly different from unity, with the largest corrections being of order 
GMq/Rq ~ 10~^ near the Sun. The amplitude of the change in 0(T) for the 1/R 
model implies a change in the metric of order ~ 1/3, which is comparable to the 
change in the metric when going from infinity to nearby the event horizon of a 
black hole. The difference is that this variation in the metric occurs in a much 
shorter scale, which must produce even larger accelerations. Something similar hap- 
pens to the models f{R) = R — p^'^'^^'^^ / R^ , for whi ch the change in 4> is of order 
^ nl{n + 2). To save those models, one could arguJ^^^^ that the density in the 
solar system is always much larger than p^, which could prevent (t){T) from reaching 
its vacuum value 00- However, this seems a very weak argument because the struc- 
ture of matter is discrete (localized wavepackets) and, therefore, one can always 
find regions in which 4i{T) takes all possible values, which should have observable 
consequences at microscopic scales. We will see later a detailed example of this in 
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Figure 2. Dependence of (f>{T) = fR{T)/fii{0) on p/pp in the 1/R model. The function smoothly 
interpolates between the two asymptotic constants 4>ao = 3/4 and <po = 1. 



Sec. 14.21 For the moment, we just conclude that to have a chance of being viable 
according to local experiments, any Palatini f{R) theory must be characterized by 
a function /(i?) such that fB,{T) is not very sensitive to density variations over the 
range of densities accessible to those experiments. From a technical point of view, 
this simply means that (f>{T) must be almost constant because then with a simple 
constant rescaling of the metric one can bring (f>{T) from (l)(T) « (/)o+corrections to 
4>{T) « 1+corrections, which turns the metric into its standard almost-Minkowskian 
form 5^1, = T^^y+corrections. In particular, if one accepts that local experiments are 
carried out in an environment with density piocai ^ Pfi, for the 1/R model <p would 
h8^= fR{T)/fnioo) = /fl(r) = 40(r)/3. 

Since for viable models we must have « 1 + with |f2(r)| <C 1, we 

can express the metric outside spherical bodies as 

gu~~l + 2U + niT) , (58) 

where U = GMq/t + K^jfr'^/Q, and ^{T) is sensitive to the sources present at 
radius r. This unusual local dependence must also be very weak, which can be 
used to impose constraints on the family of allowed Lagrangians. A detailed discus- 
sion of such constraints can be found in Ref. I53[ where the Newtonian and post- 
Newtonian limits of f{R) theories in metric and Palatini formalism was worked 
out. The results obtained in Ref. |53] using perturbative methods coincide with the 
Newtonian expansion from the exact solutions given her^, which provides an inde- 
pendent confirmation of their validity (up to Newtonian order at least) without the 



'^This new constant rescaling of the metric is equivalent to using the same system of units in 
Einstein and Jordan frames in regions where p 2> . 

'In Ref. I53l there seems to be a wrong sign in front of ^ef f ■ That error seems to be a transcription 
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complexities that the post-Newtonian analysis involves. Other approaches to the 
Newtonian limit of Palatini f{R) theories^''' ^ have reported the existence of 
Yukawa-type corrections to the usual Newtonian potential with a length scale of or- 
der 1^1/ yj Kef f. However, such terms are generically associated with propagating 
scalar degrees of freedom, which do not exist in the Palatini version of f{R) theories. 
Moreover, those Yukawa-type corrections should also appear in the post-Newtonian 
parameter 7, as it happens in the metric version of f{R) theories and scalar-tensor 
theories, which would be in conflict with experiments due to the long interaction 
range 1^1/ ^A^ff. Additionally, in the non-perturbative exact solution derived 
here and in the perturbative approach of Ref. [53] there is no trace of such Yukawa- 
type correction. A correcting term of the form Ap similar to the one denoted her^ 
by ri(T) accompanying the Newtonian potential in (|58|) was also found in Refs. l98l 
[99I but not in Ref. [97l An interesting interpretation of that term in the particular 
model f{R) = R+ X^R^ can be found in Ref. I100[ where that Palatini model was 
compared with its metric version. The metric version has a Yukawa-type correction 
of the form AV cx J (Pxp{x)e~^^~^°^^^ /\x — xq\, which for very short interac- 
tion range, A — ?► 0, leads to AV oc p(xo). This allows to see the density-dependent 
term as the limiting case of a Yukawa interaction when the interaction range is 
ultra-short. A similar interpretation is possible fo r th e behavior of the scalar curva- 
ture in the metric version of general f{R) theorieJ^Sl, We also mention that similar 
expressions for the non-perturbative eqs. (|T7|) . ([50]) . and ([?T|) have also been found 
in Refs llOll 11021 Though from the definitions in those wor ks, B{r) = e-"('') and 
i3(r)e^**-''-' = e'^^^\ one finds exact agreement with our results, an erroneous identi- 
fication in eq. (36) of Ref llOll (or equivalently eq. (19) of Ref ll02p leads to different 
conclusions. To be precise, they claimed that = e^" + e^*, whereas from their 
field equations (and ours) one finds = e"""*"^*. The discussion of the Newtonian 
limit given in Ref, 102 also assumes that f{R{T)) can be expanded around T ~ 
as /(T) fa /(O) — pdrflo ■ ■ ■, which for models such as the 1/R (see eq. (|57)) ) is not 
justified. 

Following Ref l53l and our previous discussion, one finds that the weak depen- 
dence of 0(T) on the density implies that a change relative to the value of (j) 
induced by a change Ap relative to the density p must be very small, which can be 
expressed as 



Jr dp 



< 1 o 



RfR 



1 - fn/RfRR 



< 1 , (59) 



This result together with the fact that R w n^p and fR{T) w 1 in local experiments, 



error because the perturbation equations of the Appendix have the right sign (compare them with 
the metric case). In any case, that sign is irrelevant for the conclusions of that work. 
JNote that if (j>{T) admits a perturbative expansion, then iji(T) l+Tdripy which in the Newtonian 
limit implies that Q{T) = —pOTft)- 
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can be reducetJ^ to the condition 

\RhM\ < 1 . (60) 
When appHed to models of the form f{R) = R- ^^(n+i) 

one finds that 



\RfB.- 



RR\ 



n{n + 1) 



2 \ "+1 



R 



n{n + 1) 



Plocal 



n+1 



(61) 



which is much smaher than unity ii n > —1 as long as piocai ^ Pfi ^ 10"^^ g/cm"^. 
This argument and others similar to this have been used in the literature to claim 
that this family of models are no t very constrained by local experiments, which jus- 
tified their cosmological analysiJ^^^^. However, as we pointed out above, it relies 
on the assumption that the density scale is not reachable under regular experi- 
mental conditions. By considering microscopic experiments, we will show next that 
this assumption is not correct. Local experiments, therefore, will be able to test 
the Palatini f{R) dynamics and impose tight constraints on the family of allowed 
models. 



4.2. Microscopic experiments 

Shortly after Vollick's proposal for explaining the cosmic speedup using the Palatini 
version of the 1/R model, it was claimed^^ that the model was in conflict with 
electron-electron scattering experiments. The argument goes as follows. Since the 
affine connection can be expressed in terms of the metric and the matter sources 
according to by inserting back this solution into the action, one ends up with a 
theory that has new interactions among matter fields and between the matter fields 
and the curvature. The original discussion of this problem was carried out in the 
Einstein frame representation of the theory, which apparently allows for a simpler 
interpretation of the action 

S=^J d'^xV^I^Rih) - ^p^^ + S„^[r\T)h^.,i^m] , (62) 

where (j) = <t>{T) is, in general, given by solving (p3)) . and in our particular case takes 
the form (f57|) with p replaced by T. The explicit coupling to the matter of the factor 
0(T) in Sm together with the new matter term -^r^ were inmediately interpreted 
as a clear indication that the theory should be i n conflict with particle physics ex- 
periments. This view was seriously criticizecff^ (gee also Ref. I104p because if the 
theory is analyzed in the original Jordan frame, the direct coupling of </i(T) with 
the matter action disappears and, on grounds of the Einstein equivalence principle, 
no new effects should be observed in a freely falling frame. This observation raised 
(again) a debate on the mathematical and physical consequences of working with 
different held redeflnitions and/or frames. A reanalysis of the problem-'^'^'^ then put 
forward the existence of non-perturbative couplings that prevented a consistent per- 
turbative treatment of some of the new interaction terms. The theory thus seemed 
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intractable and a definitive conclusion about its viability could not be explicitly 
worked out. 

The problem of the influence of the f{R) dynamics in microscopic systems was 
reconsiderecP^ with the focus on the possible effects that the gravitationallly- 
induced matter interactions could have on the stability of Hydrogen. Starting with 
the equation for a Dirac field in curved spac^, the corresponding non-relativistic 
Schrodinger equation for an electron in an external electromagnetic field was de- 
rived. The goal was to study if the density-dependent function (t>(T) present in the 
metric (|47p could have an effect on atoms. This question is pertinent because atoms 
are systems in which the matter density is localized around the nucleus and drops 
to zero as we move away from it. Since the density scale will be reached at some 
point, the stability and structure of atoms provides a natural laboratory to test the 
Palatini dynamics of theories sensitive to very low density scales. The fact that the 
average distance between atoms in a diluted perfect fiuid is much larger than their 
typical sizes guarantees that they can be seen as isolated systems immersed in a 
perfect vacuum, which provides a way out of the problem posed by the belief that 
local experiments are carried out within an environment whose density hides the 
presence of the modified dynamics. 

Neglecting the Newtonian potential corrections, the metric (|47| boils down to 
QfiL' = (f'^^Vf^i^- For a metric of this type, one finds the following Schrodinger-Pauli 
equation 

Er]=l [{p - el)^ - ea ■ B] + eAo\ 7] (63) 

lm + niQ J 

1 

m + mo 



ia{Vn X V) - 2ieiA ■ Vn) + - \Vn\^ + 2{Vn ■ V) + -mo)}r] 



where E is the non-relativistic energy of the electron, i] is a two-component spinor, 
771 = TO(/)~^/^, Too is a constant of order the mass of the electron to, a are the Pauli 
matrices, Aq and A are the components of the electromagnetic potential 4- vector, 
B is the external mag netic field, 17 = f ln<?!'(T), and T = —mrj^rj. If one considers 
GR, = 1, the usual Schrodinger-Pauli equation is recovered by just identifying 
Too with m. For concreteness, let us consider the 1/R modeP^. In this case, the 
function 0(T) is given in (|57p . and expressing length units in terms of the Bohr 
radius ao = 0.53 x 10^^" m, we find that p/ ~ lQ^^Pe{x), where Pe{x) — rfr] 
is the probability density of finding an electron. This expression for p/ p^ indicates 
that the electron reaches the characteristic cosmic density, p/pfi ~ 1, in regions 
where the probability density is near Pe(x) — 10^^**. In ordinary applications, one 
would say that the chance of finding an electron in such regions is negligible. In our 
case, however, that scale defines the transition between the high density p 3> 



''This analysis was carried out in the original Jordan frame to avoid the discussion generated by 
the use of the Einstein frame variables. 
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and the low density p ^ regions. In regions of high density, one finds that 
rapidly tends to a constant, (poo = 3/4, which leads to to = 2to/-\/3 and Vfi = 0. 
Identifying to — ^ \/3mo/2, Eq. (j63p reduces to the usual Schrodinger-Pauli equation 

£7]^S^^[{p~eAf-ea-B] + eAojri (64) 

In regions of low density, (j) tends to unity, Vfi — 0, and to — > to as p/p^ — > 0. As a 
result, the mass factor dividing the kinetic term is now a bit smaller (toq > to) than 
in the high density region, but the mass difference to — toq is no longer zero. This 
is the crucial point, because to — toq ~ — O.lSmo represents a deep potential well in 
the outermost parts of the atom (from r « 26ao to infinity), which has important 
consequences for its stability. In fact, if one assumes that the electron is initially 
in the ground state, the deep potential well that appears in the outer regions of 
the atom makes this state unstable and triggers a flux of probability density (via 
quantum tunneling) to those regions. Using time dependent perturbation theory, 
the half life of Hydrogen subject to this potential turns out to be 

Tff = ^ « 6 • lO^s , (65) 

which is in clear conflict with observations. From this analysis we extract several 
lessons. First, we have shown that the ultra low density scales that characterize 
models aimed at explaining the cosmic speedup can be reached in microscopic sce- 
narios. Second, we have seen that the modified gravitational dynamics of those 
models can have nontrivial effects on systems such as the Hydrogen atom. Third, 
we can use standard perturbative techniques to estimate those effects and constrain 
the models. The results obtained here for the 1/R model also provide a simple test 
to determine whether a given model is compatible with observations or not. Since 
the instability of the ground state is to a large extent due to the potential well in- 
duced by the mismatch between the values of to and toq in the low density regions, 
any f{R) model that yields a non- negligible difference 



can be automatically ruled out. In particular, for the family f{R) — i? — 

we find that fR{oo) = 1 and /i?(0) = 1 + n/{n + 2), which leads to — 



a + 2 



™o \/TTZjn ^ 1 ■ This quantity is small only if |n| <C 1, which yields A„ 



— TOon/4. For not too small n > 0, the results of Ref.^^'^^ could be directly used 
to estimate the half-life of the atom. However, for very small n, the estimation of 
the half-life should be reconsidered taking into account the contributions coming 
from the V51 terms, which were negligible transient potentials in the 1/R case. For 
negative values of n one should note that rather than a potential well, one finds a 
potential barrier, which would lower the energy of the ground state. In any case, 
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all these possibilities have observable effects and could be strongly constrained with 
current d ata. In t his sense, the variation in the energy levels of Hydrogen has been 
estimatedlfl^ilOSI 

for models in which the condition is satisfied, i.e., for models 
of the form f{R) = R{1 + e{R)) such that |e(i?)| ^ 1. Since that analysis was 
carried out in the Einstein frame and with a non-linear redefinition of the Dirac 
field, we will refrain ourselves from giving a detailed correspondence between our 
formulas and those. The strategy followed there to constrain the models consisted 
on determining how the energy of a photon released due to one transition changes 
relative to that emitted in another transition, which yields a quantity that is in- 
dependent of the electron mass and also seems to be independent of the choice of 
frame. Using data for the transitions from the initial state (n, /) — (2,0) to the 
final state {n,l) = (1,0), and from {n,l) — (8,3) to {n,l) — (2,0), the following 
constraint was found 

< 4 X 10"''" (67) 

For the family of models f{R) = R — this constraint implied that 

\n\ < 10~^^. This is the tightest constraint put so far on this family of models 
(recall that from CMB anisotropies and baryon oscillations the bound was around 
\n\ < 10~^) and puts forward the relevance of microscopic experiments for the un- 
derstanding of the dynamics of Palatini theories. 

The analysis of RefHSl also raised doubts about the applicability of the Palatini 
field equations to describe macroscopic systems. A careful analysis of such prob- 
lem has been explicitly carried out in Refs ll07[ 11081 (see also Ref ll09l for a related 
discussion). It was concluded there that at the classical level the physical masses 
and geodesies of particles, cosmology, and astrophysics in Palatini modified gravity 
theories are all indistinguishable from the results of general relativity plus a cos- 
mological constant. Part of this argument was supported by the assumption that 
isolated particles are stable and should not exhibit violations of energy and mo- 
mentum conservation. Though this could be true for certain Palatini models, the 
stability of microscopic systems can not be guaranteed in general. In particular, it 
is in clear conflict with the results presented here for the Hydrogen atom and the 
family of models f{R) = i?- Since there is a flux of probability density 
to inflnity, the energy and momentum of the system are not locally conserved. Thus, 
though GR is holographic in the sense that the equations of motions for a localized 
distribution of energy and momentum surrounded by vacuum can be derived by 
considering surface, rather than volume, integrals over curvature components, the 
instability of certain isolated systems in some Palatini f{R) models may prevent 
the interpretation of the exterior space-time as completely vacuum and, therefore. 



'To obtain this result we consider as valid the assumptions made in Ref 11071 11081 evaluate fn 
and fjiji in the vacuum value Rvac = {n + 2)^/("+^)/i^, and approximate Hq by Hq = /i^(n -|- 

2)l/(n + l)/12. 



Jr. 
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as not exactly equivalent to that of GR plus a cosmological constant. To the light of 
this, we believe that part of the conclusions of Refs ilOTl llOSI should be reconsidered. 



5. Other tests 

The previous sections have provided us with a good deal of information about the 
properties of Palatini /(i?) theories. We have contrasted the dynamics of these the- 
ories against cosmological, solar system, and laboratory data, and this has allowed 
us to impose severe constraints on the form of some families of f{R) Lagrangians. 
This exercise has been particularly useful for constraining models characterized by 
ultralow density scales. We now review other approaches followed in the literature 
to understand the viability and robustness of f{R) theories and which have raised 
interesting debates. We will begin by considering the structure of spherically sym- 
metric, static stars and then will focus on the initial value formulation of these 
theories. 



5.1. Stellar structure 

In this section we consider a problem that initially seemed to affect seriously the 
theoretical viability of all f{R) models in Palatini formalism. Using the T olma n- 
Oppenhei mer- Volkov (TOV) equations for the interior of stars in equilibriunffl^, it 
was foundmi that certain terms in those equations could blow up and form curva- 
ture singularities near the surface of spherically symmetric, static matter configura- 
tions with polytropic equation of state, p{P) = (P/A')^/^, with index 3/2 < 7 < 2. 
Since the physically interesting case 7 = 5/3 (degenerate, non-relativistic fermion 
gas) lies within this interval, this result was regarded as a serious theoretical con cern 
about the viability of Palatini /(i?) theories. The problem was soon reconsidered^^ 
and interpreted differently, claiming that it had more to do with the peculiarities 
of the equation of state used than with the own structure of Palatini f{R) theories. 
This was based on the observation that for neutron stars the tidal acceleration due 
to the surface singularity becomes equal to the Schwarzschild value of GR only at 
a distance ~ 0.3 fermi from the surface of the star, which makes unrealistic the use 
of a polytropic equation of state. However, this conclusion was also challenge 
claiming that the fluid approximation is still valid on the scales at which the tidal 
forces diverge just below the surface of a polytropic sphere in the case of th e generic 
functions f{R) considered. This debate was independently reexamined 11141 reaching 
an intermediate answer, which is the one that we present here. 

Consider a static, spherical object described by a perfect fluid, with T^,^ = 
{p-\-P)u^Ui,-\-Pg^i,. Take Schwarzschild-like coordinates, and parametrize the space- 
time line element as 

A{r) 



January 21, 2011 1:5 WSPC/INSTRUCTION FILE PA2MG'V1-1 



28 Gonzalo J. Olmo 



where A{r) = 1 — 2M{r)/r. Inserting these inputs in the field equations one 
obtains the fohowing TOV coupled equation: 



fR,r 



fR 



r ^ 
2\ Mr 
r J r 

Pr 



3 

/rA 2 
Vr 



fR,r 
fR 



A 



fR 



fR,: 



fR,i 



2r - 3M 3 fR.r 



fR fR \r{r-2M) 4 fn 



P^ 



(0) 



(0) 



(69) 
0) 
(71) 



where fu^r — drfR, and we have defined 



(p + P) 



r(r 

a{r) = {p + P) 



2M) 

fR.P 
fR 



M - 



f + n\P-p) 
fR 



31 



/3(r) = (2r) 



fR,P 
fR 



1 



3(P + P)/J^.P 
4 fR 



(72) 
(73) 
(74) 



with fn p = dpfn. One can check that these expressions recover the GR formulas 
in the limit fR — I and f — R — 2A. Given an equation of state, P = P{p) 
OT p = p{P), one can use the above formulas to compute the structure of static, 
spherically symmetric objects. To do it, one must first express the functions / and 
ffj in terms of T = — p + 3P and rewrite the radial derivatives of fp in the form 
fR,r = fR.pPr and fn^rr = fR.pPrr + fp.ppPr- One then finds 

n'^fpR 



fR,- 



fR,PP — 



RfpR - fp 
^^fpL 



(3 - PP) 



RRR 



{RfpR - fp)'- 



{i-ppf 



K'^fpR 



{RfpR - fp) 



PPP 



(75) 
(76) 



where pp = and ppp = j^p- The terms pp and ppp are the reason for the exis- 
tence of divergences near the surface of polytropes with index 3/2 < 7 < 2. This can 
be easily seen as follows. Since polytropes are characterized by p{P) = (P/A')^/^, 
one finds that pp = p/{-fP), and ppp ~ (1 ~j)p/{jP)'^, which implies that pp and 
ppp diverge as P ^ if 7 > 1 in the first case and if 7 > 1/2 in the second case. 
Therefore, if those terms do not appear in the equations multiplied by appropriate 
powers of the pressure, divergences will be unavoidable for some values of 7. Let us 
now determine the dependence on P of the various terms involved in the equations. 
From their definitions, it is easy to see that Pr""* ~ P^/'', a{r) ^ p-1+2/7^ and 
/3(r)Pi"^ - p2(-i+2/7)^ If 7 < 2, those terms decay as P ^- yielding Pr - P^/'^, 
but if 7 > 2 then they grow. The combined result for 7 > 2 gives Pr ~ p2-3/7^ 



These equations correct some transcription errors present in Ref lll4l 
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which also falls to zero near the surface. By direct computation one also finds that 
Prr, fR,pPrr, and fn^pPr are well behaved as P — > for 7 < 2. However, the term 
fR.ppPr contained in fp.rr generates a term of the form pppP^ ^ p-2+3/'i ^ which 
diverges as P ^ for 7 > 3/2 and produces the singularities reported in Refs lllll 

[ml 

Contrary to the opinions provided in Refs lllll 11131 11151 we believe that the diver- 
gences that we have found here are not due to the differential structure of Palatini 
f{R) theories. The fact that, unlike in GR, the field equations contain derivatives 
of the matter fields (via the trace T) up to second order is not the reason for the 
existence of these divergences. To see this, one should note that, as pointed out 
in Ref lll21 the divergent behavior of the term pppP^ could be cured by simply 
smoothing the behavior of pp and ppp in the outer regions of the star using a 
different equation of state. Should the divergences exist even for regular equations 
of state, then one could blame the Palatini f{R) framework for this problem but, 
in our case, the field equations are simply collaborating with the polytropic equa- 
tion of state to the development of those infinities and, therefore, they are not 
directly responsible for the existence of those divergences. One should have in mind 
that the equations of state usually respond to statistical descriptions and involve a 
number of simplifying assumptions. In fact, an accurate equation of state at labora- 
tory densities is very complicated to derive, because electrostatic interactions and 
other subtle effect s rna sk the simpler statistical properties of the idealized Fermi 
gas approximatiorP^. The polytropic equation of state should therefore not be 
used beyond its expected regime of validity. This regime, however, may depend on 
the parameters that characterize the particular Lagrangian f{R) considered. For 
instance, if one takes the model f{R) ^ R± XR^ with A of order the Planck length 
squared A ~ Zp, which defines a density scale px = (k^A)^^ ~ 2 • 10^^ g/cm^, f or a 
neutron star the divergent term begins to be non-negligible at a density of ordejni 

Ps — ^ ^ ^ ^ ^ ^ 10^^^° g/cm^, which is well below any physical density one can 
imagm But if one uses a length scale much larger than Ip, the terms responsible 
for the divergences could begin to grow in regions where the polytropic equation 
of state may still be valid-'^-'^"^ In this sense, polytropes could still be used as a 
theoretical laboratory to constrain the parameters of f{R) modelJH^. 

5.2. The Cauchy problem 

A very natural requirement of any theory of classical physics is that a sufficient set 
of initial data should be enough to determine the subsequent evolution. One then 
says that a theory possesses an initial value formulation if appropriate initial data 
(perhaps subject to constraints) can be specified such that the dynamical evolution 



"For a free electron whose wave function is spread over the entire universe, the ratio rrie/ Rfj^^^ 
is of order ~ 10~^^* g/cm^. Therefore, a simple electron would be enough to remove all the 
singularities of this type in the universe. 
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is uniquely determined. If small changes in the initial data induce small changes in 
the solution over a compact region of spacetime and if such changes do not produce 
any changes in the solution outside the ca usal f uture of this region, then the initial 
value formulation is said to be well pose GR has a well-posed initial value 
problem, which results in a stable theory with a robust causal structure. Do Pala- 
tini f{R) theories have a well-posed initial value formulation? Recent workJU^^I^ 
have concluded that, unlike their metric version, Palatini f{R) theories are in gen- 
eral neither well-formulated nor well-posed, which seems a very serious reason for 
concern. We will see next, however, that Palatini f{R) theories do admit in general 
a well-formulated initial value problem. We will also use those results to argue that 
the initial value problem is likely to be well-posed. 



5.2.1. Hamiltonian formulation 

To show that the initial value problem of Palatini f{R) theories is well formulated, 
we consider the Brans-Dicke repre sentation of Pala tini the ories and work out its 
3-1-1 Hamiltonian descriptior[12l] in the usual (from now on we use 

lower-case latin letters to represent space-time indices). Consider a foliation of the 
spacetime manifold Ai into hypersurfaces St of simultaneity characterized by a 
function T(x) ^constant, a normalized timelike co-vector oc daT normal to this 
hypersurface, and a shift vector iV" orthogonal to n° = g"''«-b. This allows us to 
construct a time flow vector — Nn'^' + iV", where A'' is known as lapse, and 
decompose the metric in the form gab = hab — naUb- Elementary, though lengthy, 
manipulations allow us to express the Lagrangian density of (|29|) as follows 



where Kab — h'^hfV is the extrinsic curvature, Da4> — h\Vb4>, 4' — f^dacj), ^"^^R 
is the Ricci scalar of the 3— metric hab, and we have used the following relations 




(77) 



+ [KabK'''' - {Ka'^ f + 2Ve J^] 

nf^Va-n" - n^Van" 



(78) 
(79) 
(80) 




(81) 
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The canonical variables of the theory are {gabT4>) = {N, N°' , hab, (/>)■ The canonical 
momenta are defined by the following expressions 



n 



N 



ss_ 

SN ' 

ss 



, n. 



^5* 
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Shab 



2^ 



'ab 



LftT"" - Kh'"') U - N^D. 
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2k2 V N(j) 



N^Dr 



(82) 
(83) 
(84) 



Like in GR, we immediately see that the momenta conjugated to N and iV" are 
constrained to vanish. On the other hand, from the combination of lih = ^afcll"^ 
and TT^, we find that 



3 + 2w;\ y/h 
N ) 2^ 



N^Dr 



(85) 



is also constrained to vanish when uj — —3/2, which is the case that interests us. 
It is now useful to rewrite the Lagrangian density £ using the definition for 11'^'' to 
eliminate the explicit dependence of Kab from it. The result is 



2^ 



+ (3 + 2w) 



2N(j) 



N ~NV{(t>) 



(86) 



Note that when w ~ —3/2, the last term in the above equation vanishes whereas it 
persists for w ^ —3/2 and can be expressed in terms of the momenta using (j85p . A 
detailed discussion of both cases can be found in Ref J2li Here we will just focus on 
the Palatini case, w = —3/2. To proceed with the construction of the Hamiltonian 
one must have in mind the above constraints and apply Dirac's alg orithnll^ for 
constrained Hamiltonian systems. 



Like in GR, we have the primary constraints Cat = Hjv(i,x) = and Ca = 
Hci(i,a:) = 0. Additionally, we have the constraint ((85|) . The Hamiltonian is con- 
structed by introducing Lagrange multiplier fields \t^{t^ x\ Xait, x), and for the 
primary constraints and performing the Legendre transform as usual with respect 
to the remaining velocities. The result is 



(87) 
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where 





= Hjv , Ca 






= Ilh- (t>'^4 
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,(90) 

(91) 



For the dynamics to be consistent, the constraints must be preserved under evolu- 
tion, which requires that Cn = {H,Cn} = and Ca = {H,Ca} = 0, where the 
poisson bracket at time t is defined as 



{A{x),B{x')} = j d^a 



5A{x) 5B{x') 5B{x') 5A{x) 



(92) 



where 11* and Qi generically represent the canonical variables. By direct evaluation, 
one finds that Cn = -SH/6N = -Hn and Ca = -SH/SN" = -Ha- We thus see 
that on consistency grounds we must impose the secondary constraints Hn = and 
Ha = 0, which together with = implies that the Hamiltonian H is constrained 
to vanish, like in GR. If matter is present, one must add the corresponding pieces 
^Hmatt/^N and dHmatt/^^"' to these constraints, which leads to 



1 5Hr. 



a 5N 



■ = 

(93) 



2DJ11 + n^Da^ 



1 5Hr. 



(94) 



where we have defined a = /i^/^/(2k^) and used the tilde to denote the tensorial 
quantities tt,/, = -Kip/a and 11"^ = U^^/a. After some lengthy algebra, one finds the 
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following evolution equations 

= N^Da<l) - A 

^4, = N 



N 



dV 



2AN + 2wDc —D^cf) + N''Da^4, - 



2N 



hah^ 



hab = 2D^aNb) + —[Tiab- -f^h ] + Mhab 



(95) 
(96) 

(97) 



jjacjjb n™"n, 

4 



20 



+ D"D^(N(j3) - h''''A{N(j>) - 2D''ND^(I) + h'^'^D^ND" 
2 a (S/iab 2 



(98) 



Using these evolution equations and the constraint (|93| . one can verify that the 
evolution of leads to 

= {if, C4 = -2a7Vy(0) + a7V0^ - - /^ab^^ (99) 



2 SN 

Since Cd, must vanish, we must impose the secondary constraint 



Shab 



.dV 



2V{^) 



1 5Hr^ 



2a SN 



1 , SHmatt „ 

Na dhab 



(100) 



Using the definitions Tab 



the fact that 



2 (5£ 



„ab Lab 1 /-la 



(i°-7V°)(t''-iV''), and 



SV-r, 



yields 



SN 



and 



(/.^ - 2y(<?!)) = 



, one can verify!^ that pOO]) 



(101) 



This equation reproduces the relation (l32t when uj = —3/2 and establishes an al- 
gebraic relation between the trace of the energy-momentum tensor of matter and 
the scalar field — 4'{T). 



5.2.2. Discussion 

From the derivations of above, we see that the dynamical variables in the Brans- 
Dicke case w — —3/2 are just {hab,^"''') (plus the [qi^p^) of the matter), because 
the evolution equations for (0, tt^), as we saw above, can be combined to establish 
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the secondary constraint (llOip . The lapse, N, and shift, N°-, manifest the difFeomor- 
phism invariance of the theory and are not dynamical variables either. It is worth 
noting that the constraint (j93|) involves up to second-order spatial derivatives of 
hab (see the term ^^^R), but only first order time derivatives of it (contained in 
the momenta 11"''). However, though that constraint contains spatial derivatives of 
(j) — (f>{T) up to second order (see the term DaDbcj)), it does not contain any time 
derivative of (/"(T) because the corresponding momentum tt^ is absent in that equa- 
tion. Something analogous occurs in the vector constraint (j94p . where we can use 
the replacement tt^ = 11^ /0 to show that no extra time derivatives of the matter 
appear in the constraints. This is a very important aspect, because it means that 
the highest order time derivative of the matter fields appearing in (|93|) and (|94|) 
is the same as in GR and coincides with the highest order present in the energy- 
momentum tensor of the matter. The evolution equations also have this property. 
A glance at (|95ll98p puts forward that the evolution equations for 0, hab, and 11°'' 
do not contain the momentum tt^, while in the equation for tt^, one can reexpress 
the term A^tt^ using tt^ = Tlh/4>- Therefore, though one can find up to second-order 
spatial derivatives of 4>{T), and hence of T, there is no trace of extra time derivatives 
acting on the matter fields. The existence of second-order spatial derivatives of 0(T') 
requires an extra degree of smoothness in the matter profiles, an aspect that is not 
necessary in GR. This extra degree of differentiability is a natural requirement if we 
attend to the f{R) formulation of the w = —3/2 theory. Since the affine connection 
is compatible with a metric tab which is conformally related with the space-time 
metric gab, the smoothness and differentiability of the conformal geometry is guar- 
anteed if the conformal factor is differentiable up to second order (to yield a smooth 
field strength, Riemann tensor, of the affine connection). 

Let us now focus on the initial value problem. It is well-known that if in GR 
one specifies initial values for N, N°-, hab and 11"'' which are consistent with the con- 
straint equations, the evolution equations uniquely determine hab and 11"'', while 
N and N"^ remain undetermined, which expresses the existing gauge freedom of 
the theory. This guarantees that t he intrin sic (coordinate-independent) geometry of 
space-time is determined uniqueljl^IE^ by an initial choice of hab and 11"'' . The 
same is true for the scalar-tensor theories considered here, thus implying that the 
initial value problem is well-formulated for all w. For the w = —3/2 case, the 
only difference with respect to GR is that one must specify an initial value for 
taking into account its corresponding constraint equation to consistently establish 
the initial data. 

Though the evolution equations presented here are not suitable to determine 
whether the initial value formulation is also well-posed, it is well-known that us- 
ing different variables and representations of the evolution and constraint equations 
one can proof the well-posedness of GR an d of g eneric Brans-Dicke theories with 
w 7^ —3/2 in both Einstein and Jordan frame a 1 One can also make special choices 
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for the lapse-shift pair and manipulate the corresponding 3 + 1 equations of GR to 
show that the conjugate variables hab and U°^^ do satisfy a hyperbolic evolution 
systenff^. One can thus exploit the resemblance between the constraint and evo- 
lution equations (|93)) . ((94)) . pOTj) and ([971)- (IMl) with those of GR to argue that the 
Cauchy problem is likely to be well-posed also for the Brans-Dicke case w = —3/2. 
Note first that in vacuum, T^i, = or Hmatt = 0, the constraint (jlOip implies that 
(/) is a constant, (/)o, which turns the constraints and into 



+ V^(0o)=O (102) 

- 2DJ{1 = . (103) 

With a simple constant rescaling of the metric, these constraints are the same as 
those of GR with a cosmological constant. Setting for consistency the Lagrange 
multiplier = 0, the evolution equations for h^b and 11"'' also recover the same 
form as those of GR with a cosmological constant. We can thus conclude that the 
Cauchy problem in vacuum is well-posed. 

When matter is present, one should add to the above equations those correspond- 
ing to the matter fields. The strategy now would be to interpret the (/)-dependent 
terms, which are functions of the trace T, as part of a new (or modified) matter 
Hamiltonian. This way, the constraint and evolution equations maintain a struc- 
ture that closely resembles that of GR except by some non-constant factors ^(r) 
that multiply or divide objects like ^'^•'i? and ri''''na6. If the matter fields satisfy 
the spatial differentiability requirements imposed by the constraint equations, the 
absence of higher-order time derivatives of the matter fields in the constraint and 
evolution equations suggests that the time evolution will be as well-posed as in GR. 
This, in fact, has been explicitly shown f or a perfect fiuid using the Einstein frame 
representation of the evolution equation j-^^^ i Obviously, since in general the well- 
posedness of the GR equations depends on the particular matter sources considered, 
the modification of the source terms induced by the existence of (/)(T)-dependent 
terms requires a model by model analysis. Therefore, though one cannot conclude 
that the Cauchy problem is well-posed for an arbitrary f{R) Palatini Lagrangian, 
we find no reasons to suspect that it is ill-posed in general. 

To close this section, we comment on recent literature that criticizes the viabil- 
ity of all Palatini f{R) theories based on a seemingly ill-formulation of the Cauchy 
problem. In Ref lll9l it was claimed that the disappearance of the d'Alambertian 
□(/) from (p2|) for the value w = —3/2 implies that the non-dynamical field </> can 
be arbitrarily assigned on a region or on the entire spacetime, provided its gradient 
satisfies a degenerate equation [Eq. (4.5) in that paper], which reduces to a con- 
straint. This fact, it was stated, would make impossible to eliminate the term 
from the evolution equations unless \3<j) = 0. This was interpreted as a no-go theo- 
rem for Palatini f{R) gravity, which would have an ill- formulated Cauchy problem 



January 21, 2011 1:5 WSPC/INSTRUCTION FILE PA2MG'V1-1 



36 Gonzalo J. Olmo 

even in vacuum. This interpretation is conceptually wrong (see also Ref ll27l in this 
respect) because the scalar field in the w — —3/2 is just a given algebraic function 
of the trace T and, therefore, is clearly specified by the local matter contenl|^. More- 
over, one should note that Eq.(4.5) of Ref lll9l is not correct. That equation should 
recover the well-known relation 2V — (pV = k^T that establishes the algebraic rela- 
tion between (p and T [the secondary constraint (|10ip ]. Using Eqs. (3.4) and (3.5) of 
Ref lllQl it is easy to check that the associated Eq. (3.10) does recover our equation 
(|10ip in the Brans-Dicke case w = —3/2 (even though this is not the result obtained 
in Ref lllQp . This indicates that the first claims against the well-posedness of the 
Cauchy problem for Palatini f{R) theories stemmed from a misleading analysis of 
erroneous equations. 

The strong conclusions of Ref lllQI were a bit relaxed in Ref ll20l (see in this sense 
Refs ll29(ll30p . where it was admitted that the Cauchy problem should be well-posed 
in vacuum and with radiation fields (for which T — and (j) ^constant). In fact, 
in Ref ll20l it was correctly noticed that in the w = —3/2 case the field (j) could be 
algebraically solved in terms of T (though their Eq. (219) is the same as Eq.(4.5) of 
Ref lllQp . It was then argued that the existence of terms of the form 0(f){T), which 
imply contributions of the form DT, would cause problems for the Cauchy problem. 
Though such terms and the possible existence of higher-order derivatives of the mat- 
ter fields are certainly a reason for concern, it was prematurely concluded that the 
Cauchy problem for Palatini f{R) theories was likely to be neither well-formulated 
nor well-posed unless the trace T were constant. These conclusions contrast with 
the findings of Ref ll21l presented here, which show that the evolution equations do 
not introduce higher-order time derivatives of the matter fields, which guarantees 
that the initial value problem is as well formulate d'^ 1 1 1 22 J -^^ 

6. Nonsingular bouncing cosmologies 

We have seen that cosmological observations and local experiments strongly con- 
strain the form of the f{R) gravity Lagrangian at low curvatures (see RefslTSl [79l 
inillTJIHaillllillllllillSllOliB cosmologlcal constraints and RefsISl US US 
11061 11071 llOSI for local experiments). Though many f{R) models have the ability to 
produce late-time cosmic acceleration and fit well the background expansion history, 
they are not in quantitative agreement with the structure and evolution of cosmic 
inhomogeneities. Additionally, we have seen that the fact that matter is concen- 
trated in discrete structures like atoms causes the modified dynamics to manifest 
also in laboratory experiments, which confirms earlier suspicions on the viability 



°The fact that the ampUtude of the scalar field when ui = —3/2 is determined algebraically by the 
local matter sources also implies that the effective Newton's constant G^ff = G/4> is only subject 
to local variations of the energy-momentum densities. In this sense, though Gf,ff does change 
over cosmic timescales due to the expansion of the universe, it is not subject to the same type of 
time evolution that affects the w ^ —3/2 Brans-Dicke theories and other dynamical scalar-tensor 
theories, which contrasts with the interpretation of Ref ll28l 
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of such models according to their corresponding Newtonian and post-Newtonian 
properties. This is a very disturbing aspect of the models with infrared corrections, 
which demands the consideration of a microscopic description of the sources and 
prevents the use of macroscopic, averaged representations of the matter. A careful 
analysis of this point put forward the existence of non-perturba^^ e ffects induced 
by the Palatini dynamics in a number of context j4» l i^5 | 21 | lll | 113 | lU6 | 

In this sense, 

it is worth noting that even though the ground state of Hydrogen can be studied 
using standard perturbative methods, the first and higher excited states do man- 
ifest non-perturbative properties Despite the fact that the modified dynamics 
is strongly suppressed in regions of high density, non-perturbative effects arise near 
the zeros of the atomic wavefunctions, where the matter density crosses the charac- 
teristic low-density scale of the theory and the gradients of the matter distribution 
become very important for the dynamics [see eq. (j25p ]. Though this certainly is an 
undesired property of infrared-corrected models, it could become a very useful tool 
for models with corrections at high curvatures. Can we construct singularity-free 
cosmological models that recover GR at low curvatures using the non-perturbative 
properties of Palatini theories? As we will see, ultraviolet-corrected Palatini models 
turn out to be very efficient at removing the big bang cosmic singularity in various 
situations of interest. In this section we will thus review recent efforts carried out 
to better understand the properties of Palatini theories in the early universe. 

6.1. Non-singular f{R) cosmologies 

Growing interest in the dynamics of the early-universe in Palatini theories has 
arisen, in part , from the observation that the effective equations of loop quantum 
cosmologyl^ll (LQC), a Hamiltonian approach to quantum gravity based on the 
quantization techniques of loop quantum gravity, could be exactly reproduced by 
a Palatini f{R) Lagrangian-*^"^^. In LQC, non-perturbative quantum gravity effects 
lead to the resolution of the big bang singularity by a quantum bounce without in- 
troducing any new degrees of freedom. Though fundamentally discrete, the theory 
admits a continuum description in terms of an effective Hamiltonian that in the 
case of a homogeneous and isotropic universe filled with a massless scalar field leads 
to the following modified Friedmann equation 



where pcrit ~ 0.41pp;ancfc- At low densities, p/pcrtt ^ 1, the background dynamics 
is the same as in GR, whereas at densities of order pcrit the non-linear new matter 
contribution forces the vanishing of and hence a cosmic bounce. This singularity 
avoidance seems to be a generic feature of loop-quantized universei^^. 
Palatini f{R) theories share with LQC the fact that the modified dynamics that 
they produce is not due to the existence of new dynamical degrees of freedom 
but rather to non-linear effects produced by the matter sources, which contrasts 




(104) 
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with other approaches to quantum gravity and to modified gravity. This similarity 
makes it tempting to put into correspondence Ea. (|104p with the corresponding f{R) 
equation 

fR {k^p + {RfR - /)/2) 



(105) 



Taking into account the trace equation ((22|) . which for a massless scalar becomes 
RJr — 2f — 2k^p and implies that p — p{R), one finds that a Palatini f{R) theory 
able to reproduce the LQC dynamics (|104l) must satisfy the differential equation 

A/r-B 



miR ~ 3/)A + RB 



(106) 



where A - ^2{RfR - 2/)(2i?, - [RJr - 2/]), B = 2^RjR{2RfR - 3/), and 
Rc = K^Pc- If one imposes the boundary condition Iutir^o Jr. — > 1 at low cur- 
vatures, and Hlqc = apai (where a represents the acceleration of the expansion 
factor) at p = p r„ the solution to this equation is unique. The solution was found 
numerically!^, though the following function can be regarded as a good approxi- 
mation to the LQC dynamics from the GR regime to the non-perturbative bouncing 
region (see FigE]) 



dR 



tanh 



103 



In 



R 



l27^, 



(107) 



It should be noted that different attempts to find effective actions for the LQC 




Figure 3. Vertical axis: df /dR ; Horizontal axis: R/Rc- Comparison of the numerical solution 
with the interpolating function I I107I I. The dashed line represents the numerical curve. 



equations have also been c onsid ered but either failecli:^ or are limited to the low- 
energy, perturbative regim^l^. 
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Though the function (|107p imphes that the LQC Lagrangian is an infinite series, 
which is a manifestation of the non-local properties of the quantum geometry, the 
fact is that one can find non-singular cosmologies of the f{R) type with a finite 
number of terms. In fact, a simple quadr atic L agrangian of the form f(R ) = R + 
R'^/Rp does exhibit non-singular solutionJI^ for certain equations of stat d^'^^l^'^^ l 
depending on the sign of Rp. To be precise, if Rp > the bounce occurs for sources 
with w ^ P/p > 1/3. If Rp < 0, then the bouncing condition is satisfied by it; < 1/3 
(see Fig|4|). This can be easily understood by having a look at the expression for the 
Hubble function in a universe filled with radiation plus a fluid with generic equation 
of state w and density p 



/ + (1 + 3w)k'^P + 2K'^Prad 



[1 + |Ai] 



(108) 



where Ai = -il + w)pdpfR/ fp = il + w)il-3w)K^pfRR/{fpiRfRR~ fp)). Due to 
the structure of Ai, one can check that vanishes when /^j — >■ 0. A more careful 
analysii^ shows that /i^ — is the only possible way to obtain a bounce with a 
Palatini f{R) theory that recovers GR at low curvatures if w is constant. In the 
case of f{R) = R + R^ /Rp, it is easy to see that /r = has a solution if 1 + 

'^Rsounce/Rp = is Satisfied for ppounce > 0, where Rpounce = (1 - iw)n'^ p Bounce, 

which leads to the cases mentioned above. It is worth noting, see Fig|31 that the 
expanding branch of the non-singular solution rapidly evolves into the solution 
corresponding to GR. The departure from the GR solution is only apparent very 
near the bounce, which is a manifestation of the non-perturbative nature of the 
solution. Note also that in GR there is a solution that represents a contracting 
branch that ends at the singularity where the expanding branch begins (this solution 
is just the time reflection of the expanding branch). The Palatini model f{R) — 
R—R^/2Rp represented here simply allows for a smooth transition from the initially 
contracting branch to the expanding one. 

Besides avoiding the development of curvature singularities, bouncing cosmolo- 
gies can solve the horizon problerrP^ which makes them interesting as a substitute 
for inflation. To be regarded as a serious candidate to explain the phenomenology of 
the early universe, these theories should provide a consistent evolution of perturba- 
tions across the bounce, which should also be compatible with the observed nearly 
scale invar iant spectrum of primordial perturbations. Investigations in this direction 
have founcPSni ^]2at /(i?) models that develop a bounce when the condition fR = 
is met turn out to exhibit singular behavior of inhomogeneous perturbations in a 
flat, dust-filled universe. However, since some terms in the perturbation equations 
blow up as fp — >■ 0, their backreaction renders the perturbative system invalid and, 
therefore, one cannot say if there is a true singularity or not. 



Furt her insight on the robustness of the bounce under perturbations was 
obtainecP^ studying the properties of f{R) theories in anisotropic spacetimes of 
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GR -Vs- f (R)= R 

a[t] 



f (R) , K>0 

- - - GR , K>0 
f (R) , K=0 

- - - GR , K=0 
f (R) , K<0 

- - - GR , K<0 



-300 -250 -200 -150 -100 -50 ' 

Figure 4. Time evolution of the expansion factor for the model f{R) = R — /2Rp and w = 
for K > 0, K = 0, and K < (solid curves from left to right). From left to right, we see that the 
universe is initially contracting, reaches a minimum, and then bounces into an expanding phase. 
The dashed lines, which are only discernible near the bounces, represent the expanding solutions 
of GR, which begin with a big bang singularity (a(t) = 0) and quickly tend to the nonsingular 
solutions. 




Bianchi-I type 

ds^ = -df + > ^ ai{t){dx' f . (109) 



3 

1=1 



If one considers these space-times under the dynamics of Palatini theories with a 
generic perfect fluid, one can derive a number of useful analytical expressions. In 
particular, one finds that the expansion 6 — Hi and the shear cr^ = (iJ^ — |) 
(a measure of the degree of anisotropy) are given by 



= -F2 o . (Ill) 



where the constants Cij = ~Cji set the amount and distribution of anisotropy and 
satisfy the constraint C12 + C23 + C31 = 0. In the isotropic case, Cy = 0, one 
has cr^ = and 6*^ = 977^, with given by Eq. (|108p . Now, since homogeneous 
and isotropic bouncing universes require the condition /^^ = at the bounce, a 
glance at ()111|) indicates that the shear diverges as ~ !//!;■ This shows that, re- 
gardless of how small the anisotropics are initially, any isotropic f{R) bouncing 
model will develop divergences when anisotropics are present. It is worth not ing 
that even though diverges at = 0, the expansion and its time derivative^ 
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are smooth and finite functions at that point if the density and curvature are 
finite. However, one can check by direct calculation that the Kretschman scalar 
Rt.uapR'"""' = 4.{J2^{H^ + Hf)^ + + Hi + Hi) diverges at least as 

~ which is a clear geometrical pathology and signals the presence of a phys- 

ical singularity. The problems when /r v anishes seem to be generic in anisotropic 
models of modified theories of gravity^^. 



6.2. Nonsingular cosmologies beyond f{R) 

The previous section provides reasons to believe that Palatini f{R) models are 
not able to produce a fully satisfactory and singularity-free alternative to GR in 
idealized universes filled with a single perfect fluid with constant equation of stat^. 
Though the homogeneous and isotropic case greatly improves the situation with 
respect to GR, the existence of divergences when anisotropics and inhomogeneities 
are present spoil the hopes deposited on this kind of Lagrangians. To the light of 
these results, new Palatini theories were explorecP^to determine if the introduction 
of new elements in the gravitational action could avoid the problems that appear 
in the f{R) models. This led to the study of isotropic and anisotropic cosmologies 
of some simple generalization of the f{R) family in which the Lagrangian takes the 
form f{R, Q), with Q — Rfj^^R'^'^. Using the particular Lagrangian 

/(i?,i?^.i?'''') = i? + a|^ + %^ , (112) 
-ftp -ftp 

where Rp ^ Ip^ is the Planck curvature, it was found that completely regular 
bouncing solutions exist for both isotropic and anisotropic homogeneous cosmolo- 
gies filled with a perfect fluid. In particular, one flnds that for a < the interval 
0<w<l/3is always included in the family of bouncing solutions, which contains 
the dust and radiation cases. For a > 0, the fluids yielding a non-singular evolu- 
tion are restricted to w > ^^^^ , which implies that the radiation case it; = 1/3 is 
always nonsingular. For a detailed discussion and classification of the non-singular 
solutions depending on the value of the parameter a and the equation of state w, 
see RefEll 

The field equations that follow from the Lagrangian (|112p when -R^^ is assumed 
symmetric!^ were derived in ReflHS! (see also Refs ll07l I143P and take the form 

fR.R^j.v — ^Qt^u + "^fgRf^aR^iy = K^T^jy , (113) 
[V^ ifRg^" + 2fQRn] = (114) 



PThe consideration of se vera l fluids, fluids with varying equation of state, or fluids with anisotropic 
stresses, see for instancJHI] could affect the dynamics providing new bouncing mechanisms and 
preventing the extension of this conclusion to such more realistic cases. 
''See Ref ll42l for the case when this condition is relaxed. 
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where fu = ObJ and Jq = dqf. The connection equation (|114l) can be solved 
in general introducing an auxiliary metric hap such that (|114p takes the form 



0, which implies that F^^^^ can be written as the Levi-Civita 



connection of hnp. When the matter sources are represented by a perfect fluid, 



T^u = (p + P)Uf^u^ + Pg^,^, one can show that h^^ and h'^'^ are given hy 



56_. 



1 
n 



A. 



fj.i' 



Ai- 

A2 



A. 



where 



f2 = [Ai(Ai-A2) 



Ai 



1/2 



f2 
J R 



Ao 



A± - + 



(115) 
(116) 

(117) 
(118) 



It is worth noting that (jll5l) implies a disformal relation between the metrics g^^i, 
and hu„. A relati on of this form between two metrics naturally arises in Bekenstein's 



relativistic theorjU44| Qf MOND and in previous versions of it. In the MOND theory, 
the vector is an independent dynamical vector field and the functions in front 
of it and in front of (7^,^ depend on another dynamical scalar field. In the theory 
described here, on the contrary, the metric tensor is the only dynamical field of 
the gravitational sector. Note also that a Palatini-like version of MOND has been 
recently proposed by Milgro 

In terms of /i^ii/ and the above definitions, the metric field equation ()113p takes the 
following form 



Ai 



{f + 2K^P) AlK^p + P) 



2n 



Ai - A2 



(119) 



In this expression, the functions /, Ai, and A2 are functions of the density p and 
pressure P. In particular, for our quadratic model one finds that R — K^{p — 3P), 
like in GR, and Q — Q{p, P) is given by 



Q 

2Rp 



/ 

2 



Rp 
32 



fn.] - 



An^{p + P) 
(120) 



where / — R+aR^ / Rp, and the minus sign in front of the square root has been cho- 
sen to recover the correct limit at low curvatures. In a universe filled with radiation, 
for which R = 0, the function Q boils down tJ^ 



1-^-, 1 

3Rp V 



16k?p 
3Rp 



(121) 
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This expression recovers the GR value at low curvatures, Q « 4(k^p)^/3 + 
32{k^p)^ /9Rp + . . . but reaches a maximum Qmax = 3Rp/16 at n^Pmax — 3i?p/16, 
where the squared root of (|120p vanishes. At pmax the shear also takes its maximum 
allowed value, namely, 



: A/3/16i?p ^(Ci2 + + C|i), which is always finite, 
and the expansion vanishes producing a cosmic bounce regardless of the amount 
of anisotropy (see FiglS]). The model j[ll2l) . therefore, avoids the well-known prob- 
lems of anisotropic universes in where anisotropics grow faster than the 
energy density during the contraction phase leading to a singularity that can only 
be avoided by sources with w > 1. 



Q 

f(R,Q)= R+a— + ^ 

Kp Kp 




0.05 



0.10 



0.15 



Figure 5. Evolution of the expansion as a function of K^p/Rp in radiation universes witii low 



anisotropy, which is controlled by the combination = 
corresponds to the isotropic flat case, 9^ = 9H^. 



' '^23 + *^31 • 



The case with C 



: 



The evolution of inhomogeneities in the quadratic model discussed here was con- 
sidered in Ref ll46l though the approximations used there to solve for the connection 
equation did not allow to see the existence of bouncing solutions. For this reason, in 
this case one cannot make any statement regarding the evolution of inhomogeneities 
across the bounce. The cosmology of f{R) and f{R^vR^^) theories was also con- 
sidered in some detail in Ref ll47l The possibility of having a standard cosmological 
evoluti on in /(i?, Q) models with a large cosmological constant has been considered 
recentljffSHI 

It should be noted that the choice of a symmetric Ricci tensor in the analysis of 
/(i?, Q) bouncing cosmologies presented above is not arbitrary. As shown in Ref ll42[ 
the antisymmetric part of the Ricci tensor introduces new dynamical degrees of free- 
dom in the form of a massive vector field (see also Ref ll49l for a related result). If 
one looks for a framework suitable for the description of the effective dynamics of 
LQC (including anisotropics) and, more generally, of other theories of quantum ge- 
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ometry not involving new degrees of freedom, it seems natural to impose constraints 
on the spectrum of possible Lagrangians to avoid new propagating fields. In this 
sense, we note that the /(-R, Q) theories discussed here are able to reproducJlSO] 
other aspects of the expected phenomenology of quantum gravity at the Planck 
scale. In particular, without imposing any a priori phenomenological structure, the 
quadratic Palatini model (I112p predicts an energy- density dependence of the metric 
components that closely ma tches the structure conjec tured in models of Doubly (or 
Deformed) Special Relativitjff^ and Rainbow Gravitjff^. This confirms that Pala- 
tini theories represent a new and powerful framework to address different aspects 
of quantum gravity phenomenology. 



7. Summary and Conclusions 

From the number of works that have been discussed in this review, it seems fair to 
say that the Palatini approach to modified gravity has experienced a recent period 
of accelerated expansion motivated by theoretical and observational advances in 
cosmology. The possibility of explaining the cosmic speedup problem in geometrical 
terms boosted the interest in all sorts of modified theories of gravity with special 
emphasis in the f{R) family. Palatini f{R) theories appeared at first as an exotic 
alternative to the more familiar metric formulation of those theorie s. T hey had the 
advantage of naturally producing an effective cosmological const ant!^ , of avoiding 
certain dynamical instabilities present in their metric formulatiorl^ESl^ and of yield- 
ing second-order evolution equations. However, the first models chosen to attack 
the cosmic acceleration problem (see Secl3]) had the undesired feature of requiring 
a microscopic description of the matter sources (see Sec l4.2l) . The analysis of the 
dynamics of those models in the microscopic world put forward the existence of 
non-perturbative effects which seemed to be in clear conflict with our understand- 
ing of the physics at small scales. To overcome the technical difficulties posed by 
this situation, different directions were followed to test the viability of various fam- 
ilies of f{R) models. This motivated the analysis of the weak field l imit^— ST QS> 
the subtleties in the description of averag ed distributions of rtiatteJ^ ^ j -*^^^ ! -*^^^ ! ^^^ ! 
the stab ility and stru cture of stellar objectil^Qt l^^ l ^^^ l ^^^ l ^^'^ l ^^^ l ^^'^ l, ^j^g Cauchy 
probleml^3l2IIl26|^ and other issues that complemented the continuous investiga- 
tion of the cosmological dynamics of these theories. All these different approaches 
have raised interesting and healthy debates that have shed light on the many scenar- 
ios in which the gravitational dynamics of Palatini theories may have an infiuence. 
From those debates it follows that the background expansion history of viable f{R) 
models is currently (statistically) indistinguishable from that of the standard ACDM 
model (GR with cold dark matter and a cosmological constant), that laboratory and 
solar system tests can efficiently put constraints on model parameters, that stellar 
structure can also be used to set some constraints on f{R) models, that the Cauchy 
problem is well- formulated and well-posed in many situations of interest, and that 
Palatini theories are a powerful new tool to address different aspects of quantum 
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gravity phenomenology. 

The observation that Palatini theories can be used t o des cribe the effective ge- 
ometry of space-times with a discrete quantum structurJI^ provides solid reasons 
to explore the properties of f{R) and more general theories in the early universe 
and in scenarios involving strong gravitational fields and very high energy densities. 
In this sense, we note that simp le extensions of the f{R) family that include Ricci 
squared termPEHHHSESl 

present a much richer phenomenology than f{R) 
models. On the other hand, the fact that the dark matter problem in galaxies can 
be addressed from a new clasJ^ of Palatini theories suggests that new approaches to 
the dark matter and dark energy problems beyond the f{R) family are possible. The 
exploration of the field equations, cosmology, black hole formation, stellar structure, 
galactic dynamics, ... of new and more general Palatini theories will surely yield in- 
teresting new results with potential applications to quantum gravity, the late-time 
cosmology, and astrophysics. We hope that this review helps active researchers in 
this field and encourages newcomers to continue the exploration of the Palatini ap- 
proach to modified gravity to address and solve some of the important problems 
that cosmology faces nowadays. 
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